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Abstract In earlier work, Chekhov and Fock have given a quantization of Te- 
ichmiiller space as a Poisson manifold, and the current paper first surveys this 
material adding further mathematical and other detail, including the underlying 
geometric work by Penner on classical Teichmiiller theory. In particular, the ear- 
lier quantum ordering solution is found to essentially agree with an "improved" 
operator ordering given by serially traversing general edge-paths on a graph in the 
underlying surface. Now, insofar as Thurston's sphere of projectivized foliations 
of compact support provides a useful compactification for Teichmiiller space in the 
classical case, it is natural to consider corresponding limits of appropriate operators 
to provide a framework for studying degenerations of quantum hyperbolic struc- 
tures. After surveying the required background material on Thurston theory and 
"train tracks", the current paper continues to give a quantization of Thurston's 
boundary in the special case of the once-punctured torus, where there are already 
substantial analytical and combinatorial challenges. Indeed, an operatorial version 
of continued fractions as well as the improved quantum ordering are required to 
prove existence of these limits. Since Thurston's boundary for the once-punctured 
torus is a topological circle, the main new result may be regarded as a quantization 
of this circle. There is a discussion of quantizing Thurston's boundary spheres for 
higher genus surfaces in closing remarks. 



1 Introduction 

One manifestation of Teichmiiller space in contemporary mathematical physics is as the 
Hilbert space and the algebra of observables for three-dimensional (3D) quantum gravity 
since E. Verlinde and H. Verlinde have argued that the classical phase space of Einstein 
gravity in a 3D manifold is the Teichmiiller space of its boundary. (Analogously, the 
classical phase space for 3D Chern-Simons theory is the moduli space of flat connections 
on the boundary; this theory was quantized in [3 IS]-) Teichmiiller space possesses its 
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canonical (Weil-Petersson) Poisson structure, whose symmetry group is the mapping 
class group of orientation-preserving homeomorphisms modulo isotopy. The algebra of 
observables is the collection of geodesic length functions of geodesic representatives of 
homotopy classes of essential closed curves together with its natural mapping class group 
action. 

Given this Poisson structure, one can turn to the problem of quantizing it, thereby 
obtaining a variant of the quantum 3D gravity description. According to the correspon- 
dence principle: (1) the algebra of observables of the corresponding quantum theory is 
the noncommutative deformation of the *-algebra of functions on it governed by the Pois- 
son structure; (2) the Hilbert spaces of the theory are the representation spaces of these 
*-algebras; and (3) the symmetry group acts on the algebra of observables by automor- 
phisms. Under the assumption that the quantization of a Poisson manifold exists and is 
unique, to solve this problem it suffices to construct a family of *-algebras, which depend 
on the quantization parameter h, and an action of the mapping class group on this family 
by outer automorphisms, and to show that the algebra and the action thus constructed 
reproduces the classical algebra, the classical action, and the classical Poisson structure 
in the limit h — *■ 0. This program has been successfully performed in -4] and |5,, and 
we describe it in Sees. 2 and 3 of this paper. Actually, the problem that was solved in 
these papers differs slightly from the original formulation because the methods of [4] are 
suited for describing only open surfaces (surfaces with nonempty boundary, components 
of which can, however, reduce to a puncture). The corresponding Teichmiiller space has a 
degenerate Weil-Petersson Poisson structure, while the mapping class group acts as sym- 
metry group. We describe the deformation quantization of the corresponding Teichmiiller 
space, the action of the mapping class group by outer automorphisms, the representa- 
tions of the observable algebra, and the induced action of the mapping class group on the 
representation space following [4J. 

As in |1 3 , the representation space of the observable algebra can also be interpreted 
as the space of conformal blocks of the Liouville conformal field theory. This program is 
under current development (see [6 ). Our construction can be therefore interpreted as the 
construction of the conformal block spaces and the mapping class group actions for this 
CFT. 

The key point of the quantization procedure is constructing quantum mapping class 
group transformations that define in a consistent way the morphisms between quantum 
*-algebras simultaneously preserving the quantum geodesic algebra. The main mathemat- 
ical ingredient of the construction is a version of the quantum dilogarithm by L. D. Fad- 
deev |T3] • We interpret the corresponding five-term relation as the only nontrivial relation 
in a certain groupoid that has the mapping class group as maximal subgroup. A similar 
construction has been made independently and simultaneously by R. M. Kashaev [5]. The 
key difference between these two constructions lies in the dimensions of the Poisson leaves 
of the two theories: given a graph with v three-valent vertices, e edges, and / = s faces 
suitably embedded in a surface F* of genus g with s ideal boundary components, the 
genus is given by Euler's formula v — e + f = 2 — 2g. In Kashaev's approach, there are 2e 
variables and v + / central elements (at each vertex and at each face), so the Poisson leaf 
dimension is 2e — v — f = 8g — 8 + 3s, while in our approach, there are e variables and 
/ central elements, so the Poisson dimension e — f = Qg — 6 + 2s exactly coincides with 



the dimension of the Teichmiiller space of Riemann surfaces of genus g with s punctures. 
The approach of j3] is thus appropriate for describing 2D topological theories while that 
of is suitable for describing Liouville field theory as a lattice theory (for instance, the 
Liouville field central charge can be calculated, see [T3]). 

One of the mathematical tools employed in the quantization jl] is the decorated Te- 
ichmiiller theory [7j, and the relevant aspects are briefly reviewed in Section 2. In effect 
in the classical case, to each edge of a trivalent "fatgraph" F (i.e., a graph plus a cyclic 
ordering of the half-edges about each vertex) embedded as a "spine" (i.e., a deformation 
retract of the surface) F = Fjj is assigned a number Z a , where a here and below indexes 
the edges of T. The tuple (Z a ) gives global coordinates on an appropriate Teichmiiller 
space T H = T H (F) of F (as first studied by Thurston and later by Fock), where any or 
all of the "punctures" of Fjj are permitted to be uniformized instead as circular boundary 
components (and the subscript H stands for "holes"); the details are given in Section2.1. 
The Weil-Petersson Kahler two-form, which is known in the (Z a ) coordinates, pulls 
back to a degenerate two-form, i.e., to a degenerate Poisson structure on Th- Further- 
more, the action of the mapping class group MC = MC(F) of F on which is again 
known [Hj in the (Z a ) coordinates as well as combinatorially ^Uj, leaves invariant this 
Poisson structure, and this is the Poisson manifold Th with MC-action which has been 
quantized. 

More explicitly still, for each edge a, we may associate a pair of Mobius transformations 
R Za ,L Za depending upon Z a with the following property. For any homotopy class 7 of 
geodesic in F with corresponding closed edge-path P on T, consider the serial product P 7 
of operators Rz a ,Lz a taken in order as one traverses P in some orientation from some 
starting point, where one inserts the former (or latter, respectively) if immediately after 
traversing edge a, then P turns right (or left) in T; thus, the combinatorial geometry 
of serially traversing edges of F as dictated by P determines an ordered product P 7 of 
matrices depending upon (Z a ), and the length Z 7 of the geodesic representative of 7 for 
the point of T H determined by (Z a ) is given by G 7 = 2 cosh (Z 7 /2) = |tr P 7 |. It is the 
Poisson algebra of these geodesic functions G 7 , the algebra of observables, which has been 
quantized. 

In the quantum case (after passing to a suitable subspace on which the Poisson struc- 
ture is non-degenerate), standard techniques of deformation quantization produce the 
appropriate Hilbert space TC of the quantum theory, as well as pairs of operators Rz a , Lz a 
on Ti 2 for each edge a of F. Again, to the homotopy class of a geodesic 7 in F or its 
corresponding closed edge-path P on T, we may assign the ordered product P 7 of these 
operators as dictated by the combinatorial geometry of P in F, whose trace G 7 = tr P 7 
is the "quantum geodesic operator." The main point is to prove the invariance under the 
action of MC, which is intimately connected with functional properties of the quantum 
dilogarithm as was mentioned before. 

In jl] was proved the existence and uniqueness of an appropriate "proper quantum 
ordering" of operators that enjoy MC-invariance as well as satisfy the standard physical 
requirements. In Section 3.5, we observe that the natural operatorial ordering given by the 
combinatorial geometry of edge-paths in F can be used to derive this physically correct 
quantum ordering. The improved ordering is required in the subsequent quantization 
(discussed below). 



In order to explain the further new results in this paper, we must recall aspects of 
Thurston's seminal work on surface geometry, topology, and dynamics from the 1970- 
1980's, which is surveyed in Section 4. Very briefly, Thurston introduced the space VJ 7 o = 
VJ-"o(F) of "projective measured foliations of compact support in F" as a boundary for 
the Teichmiiller space T = T(F), where "Thurston's compactification" T = TUVJ r o is a 
closed ball with boundary sphere VJ-'o, where the action of MC on T extends continuously 
to the natural action on T. Furthermore, the sphere VJ-'o contains the set of all homotopy 
classes of geodesies in F as a dense subset, i.e., VTo is an appropriate completion of 
this set. (Unfortunately, the action of MC on VJ-'o has dense orbits, so the beautiful 
structure of Thurston's compactification does not descend in any reasonable way, with 
the current state of understanding, to a useful structure on the level of Riemann's moduli 
space.) Thurston also devised an elegant graphical formalism for understanding VJ-'o 
using "train tracks", which are graphs embedded in F with the further structure of a 
"branched one-submanifold" (cf. Section 4.2). In effect, a maximal train track gives a 
chart on the sphere VJ-'o, and furthermore, the combinatorial expression for inclusion of 
charts effectively captures the dynamics of the action of diffeomorphisms of F (cf. Section 

4.4) . 

Thus, Thurston theory arises as a natural tool to understand degenerations in Te- 
ichmiiller space and dynamics on F, and in light of remarks above, its quantization should 
provide a natural tool for studying degenerations of 3D gravity or Liouville conformal field 
theory. Many aspects of the survey of classical Thurston surface theory (in Section 4) are 
required for our subsequent quantization. 

In any case, a mathematically natural problem armed with j3] is to "quantize Thurston 
theory" : assign operators if 7 on TC to each homotopy class 7 of geodesic in such a way that 
as 7 converges in VJ-'o to a projectivized measured foliation [J 7 , ft], then the corresponding 
operators H 1 converge in an appropriate sense to a well-defined operator on TC. 

Upon choosing a spine r C F (i.e., an embedding, up to homotopy, of a graph whose 
inclusion is a homotopy equivalence), any homotopy class of curve 7 may be essentially 
uniquely realized as an edge-path on Y. Define the "graph length" g.l.(7) to be the total 
number of edges of Y traversed by this edge-path counted with multiplicities. 

In the special case of the torus Fl, we have succeeded here (in Section 5) in showing 
that the ratio of operators 

H,= p.l.( 7 )/g.l.(7) 

converge weakly to a well-defined operator as 7 converges in VJ-q{FI), where the "proper 
length" is defined by 

p.l.( 7 )= lim - tr log 2T n (ip 7 ), 

n— >oo n Z 

with T n the Chebyshev polynomials. In fact in the classical case, p. 1.(7) agrees with half 
the length of the geodesic homotopic to 7 in the Poincare metric, and in the quantum 
case, there is an appropriate operatorial interpretation, both of which are described in 
Section 5.1. In particular, for several spines whose corresponding charts cover the circle 
PJFo(F 1 1 ), the analysis involves rather intricate estimates. This leads to a natural oper- 
atorial quantization (in Section 5.3.4) of the standard simple continued fractions, which 
are intimately connected with Thurston theory on Fl (as we shall describe in Section 

4.5) . To complete the basic theory on the torus, one would like an intrinsic operatorial 



description of the circle of unbounded operators we have constructed, as well as an an ex- 
plication of the mapping class group action on it, viz., Thurston's classification of surface 
automorphisms. 

This paper is organized as follows: Section 2 covers classical Teichmiiller theory and 
Section 3 the Chekhov-Fock quantization; Section 4 surveys classical Thurston theory of 
surfaces, and Section 5 gives our quantization of Thurston's boundary for the punctured 
torus. Excerpts of Sections 2 and 3 are derived from an earlier manuscript of Chekhov- 
Fock, and we strive to include further mathematical detail. Section 4 surveys aspects 
of a large literature on Thurston theory and train tracks, explicitly covering only what 
is required in Section 5. Section 5 should be regarded as work in physics in the sense 
that some of the formal calculations depend upon manipulations of asymptotic spectral 
expansions for which there may be remaining mathematical issues. Closing remarks in 
Section 6 discuss the natural extension of these results to more complicated surfaces as 
well as other related work. Appendix A includes a novel proof of the required convergence 
in the classical case for any surface (in a sense, a new proof of the existence of Thurston's 
compactification), which may yet be useful in the quantum case. Appendix B contains 
an analysis of the Casimir operators in the Poisson algebra in the Z-variables and the 
appropriate diagonalization of Poisson structure. 
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2 Classical Teichmiiller spaces 

To begin, we shall briefly recall the two related roles played by graphs in Teichmiiller 
theory as both aspects will be required here. 

A fatgraph or ribbon graph is a graph T together with a cyclic ordering on the half-edges 
incident on each vertex, and we canonically associate to T a surface F(T) with boundary 
obtained by "fattening each edge of the graph into a band" in the natural way; we shall 
tacitly require all vertices to have valence at least three unless stated otherwise, and we 
shall call a fatgraph cubic if each vertex has valence three. To each homotopy class of 
homotopy equivalence l : T — ► F for some surface F, where i respects the orientation, 
there is a corresponding cell in T g s x R> as explained in ^T] in the hyperbolic setting 
[7j and in the conformal setting [21]. Thus, a homotopy class of i : T — > F is the name 
of a cell in the canonical cell decomposition. We shall sometimes suppress the mapping 
i : T — > F and refer to T itself as a spine of F, where T is identified with t(T) C F. 
The cell decomposition is invariant under the action of the mapping class group (induced 
by post-composition of i with homeomorphisms) , and this has been an effective tool for 
studying Riemann's moduli space; for instance, we shall recall here the corresponding 
presentation of the mapping class groups. 

The second role of fatgraphs is exclusive to the hyperbolic setting, namely, fatgraphs 
provide a kind of "basis" for geometrically natural global parameterizations of Teichmiiller 



space. Specifically, fix a homotopy class t : V — > F as above, where we now demand that 
T is cubic, and let E = E(T) denote the set of edges of T. In several different contexts, 
one can naturally identify IR> with a suitable modification of an appropriate Teichmiiller 
space of F. 

For instance, for punctured surfaces, recall [7] that the lambda length of a pair of 
horocycles is v2e*, where 5 is the signed hyperbolic distance between the horocycles. 
Lambda lengths give a global real-analytic parametrization of the decorated Teichmiiller 
space as the trivial bundle T g = T g x IRf^ over Teichmiiller space, where the fiber over 
a point is the space of all s-tuples of horocycles in the surface, one horocycle about each 
puncture (parameterized by hyperbolic length). 

For another example, Thurston's shear coordinates JH], JE] give global parameters not 
only on Teichmiiller space (cf. Section 2.1.1) but also on the related space of measured 
foliations (cf. Section rO?j) . 

On the level of Teichmiiller space, the two global coordinate systems (lambda lengths 
and shear coordinates) are closely related, and we choose to give the exposition here 
principally in shear coordinates with the parallel lambda length discussion relegated to a 
series of ongoing remarks. On the other hand, certain proofs of identities involving shear 
coordinates are easy calculations in lambda lengths. 

2.1 Graph description of Teichmiiller spaces 
2.1.1 Global coordinates on Teichmiiller space 

In addition to the Teichmiiller space T g and decorated Teichmiiller space T g s , we shall also 
require the following modification. Given an open Riemann surface F of finite topological 
type, a neighborhood of an ideal boundary component is either an annulus or a punctured 
disk; in the former case, the ideal boundary component will be called a "true" boundary 
component and in the latter will be called a "puncture." We shall study the latter as a 
degeneration of the former with an elaboration 



of the usual Teichmiiller space, where Horn' denotes the space of all discrete faithful 
representations with no elliptic elements, i.e., |tr (p(7))| > 2 for all 7 G n 1 (F) for any 
representation p : 7Ti(F) — > PSL 2 {R). 

Assume that 1 : T — > F is a homotopy class of homotopy equivalence and 7 G 7i"i (F) 
is conjugate in iri(F) to the boundary of a regular neighborhood of an ideal boundary 
component of F(T). Thus, the ideal boundary component is a puncture if and only 
if p{^j) is a parabolic transformation, i.e., |tr (p(t)) I = 2. For any 7 G t^i{F) with 
I tr (p(j))\ > 2, the underlying free homotopy class of unbased curves contains a unique 
hyperbolic geodesic whose length Z 7 is given by G 7 = 2cosh(Z 7 ) = |tr p(j)\, where G 1 
is called the geodesic operator and is constant on the conjugacy class of 7 G ^i{F). 
Furthermore, by definition, 

T g = {[p] G Teich(F) : V7 G tti{F) freely homotopic into the boundary, 




we have |tr (p(7))| = 2} 



C Teich(F) 



Finally, define the space T H (F) to be the 2 s -fold cover of Teich(F) branched over T g s , 
where the fiber is given by the set of all orientations on the boundary components of F. 

Theorem 2.1 Fix any spine Y C F , where Y is a cubic fatgraph. Then there is a real- 
analytic homeomorphsim R E ( r ) — > T H (F). The hyperbolic length Z 7 of a true boundary 
component 7 is given by Z 7 = \ where the sum is over the set of all edges traversed 

by 7 counted with multiplicity. Furthermore, ^ Zi = if and only if the corresponding 
ideal boundary component is a puncture, so T g s C Th(F) is determined by s independent 
linear constraints. 

The theorem is due to Thurston with a systematic study by Fock. We shall not give a 
proof here (though there is not a complete proof in the literature so far as we know), but 
we shall at least give the construction that defines the homeomorphism IR £(T ) — > T H (F). 

The basic idea is to associate to each edge of Y an appropriate cross ratio. To set 
this up, consider the topological surface F + C F obtained by adjoining a punctured 
disk to each true boundary component of F. The fatgraph l{Y) C F C F + is thus also 
a spine of F + , and its Poincare dual in F + is an ideal triangulation A of F + (i.e., a 
decomposition into triangles with vertices among the punctures). In the universal cover 
of F + , each arc of A thus separates two complementary triangles which combine to give 
a topological quadrilateral, and the basic idea is to associate to each edge the cross-ratio 
of this quadrilateral. 

To make this precise and describe the homeomorphism in the theorem, let a = 
\,...,E = E(Y) index the edges of Y, and let (Z a ) denote a point of R E . We asso- 
ciate the Mobius transformation 

to the edge a. To explicate this definition, consider an ideal quadrilateral in the hyperbolic 
plane triangulated by a diagonal into two ideal triangles T l5 T 2 . We may conjugate in 
PSL 2 (M.) to arrange that the vertices of Ti are 0, —1, 00 and the vertices of T 2 are 0, —1, t, 
where < t < 00, and an appropriate cross ratio of the original quadrilateral is t. Setting 
Z a = log t in the formula above, Xz a is the Mobius transformation interchanging 0, 00 
and sending —1 to t, i.e., sending 7\ to T 2 . Notice that X% is the identity in PSL 2 ( 
so X Za also sends T 2 to 7\. 

We also introduce the "right" and "left" turn matrices 

and define the corresponding operators Rz and L z , 

-Z/2 _ Z/2 

Rz = RX z = ( q ^ ]. (2.3) 



Lz = LXz 
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(2.4) 



Consider a closed oriented edge-path P in T, where we assume that P never consecu- 
tively traverses an oriented edge followed by its reverse, i.e., there is no "turning back". 
Choosing also an initial base point on P, we may imagine the corresponding curve serially 
traversing the oriented edges of T with coordinates Z\,...,Z n turning left or right from 
Zi to Z i+ i, for i = 1, . . . , n (with the indices mod n so that Z n+ i = Z\). Assign to P the 
corresponding composition 



where the matrices Lz i or Rz i are inserted depending on which turn — left or right — the 
path takes at the corresponding stage. 

Fixing any base point, the assignment P i— > Pz 1 ,...,p n £ -PSX 2 (1R) gives rise to a rep- 
resentation p E Th{F), and this defines the required map IR 5 ^ — > Tjj{F). Furthermore 
summarizing standard formulas and facts mentioned above, we have 

Proposition 2.2 There is a one-to-one correspondence between the set of conjugacy 
classes of elements of tti(F) and free homotopy classes of closed oriented geodesies in F. 
For any spine of F , each free homotopy class is uniquely represented by a cyclically defined 
closed edge-path P with no turning back, and the length of 7 is determined by 



By construction if P corresponds to a boundary component 7 of F(T), then the associ- 
ated matrix has the form R Zl Rz 2 ■ ■ ■ Rz n , or L Zl L Z2 ■ ■ ■ L Zn depending on the orientation. 
In this case, because all of the matrices R x (L x ) are upper (lower) triangular, formula 
(f2~H gives 



where the sign of this sum gives the orientation of the boundary component, and / 7 = 
corresponds to a puncture. This proves the assertions about boundary lengths. 

The Z-coordinates (i.e., log cross ratios) are called (Thurston) shear coordinates |16j.[TH] 
and can alternatively be defined by dropping perpendiculars from each of the two oppo- 
site vertices to the diagonal a of a quadrilateral, and measuring the signed hyperbolic 
distance Z a along a between these two projections. 

Assume that there is an enumeration of the edges of T and that edge a has distinct 
endpoints. Given a spine T of F, we may produce another spine T a of F by contracting 
and expanding edge a of T, the edge labelled Z in Figure 1, to produce r a as in the figure; 
the fattening and embedding of T a in F is determined from that of V in the natural way. 
Furthermore, an enumeration of the edges of V induces an enumeration of the edges of 
r a in the natural way, where the vertical edge labelled Z in Figure 1 corresponds to the 
horizontal edge labelled —Z. We say that T a arises from T by a Whitehead move along 
edge a. We shall also write T a p = (r^)^, for any two indices a,j3 of edges, to denote the 



Pzl..z, 



Lz n L Zn . 1 Rz n - 2 ■ ■ ■ RzzLzn 



(2.5) 



G 7 = 2cosh(/ 7 /2) = \ trP Zl ...z, 



(2.6) 




(2.7) 



result of first performing a move along a and then along /3; in particular, T aa = T for any 
index a. 

A B ^ 

A + cf)(Z) B-(f>(-Z) 



D-<f>(-Z) C + (j>(Z) 

D 

Figure 1- Whitehead move on shear coordinates 



Proposition 2.3 j3] Setting 4>(Z) = log(e z + 1) and adopting the notation of Figure 1 
for shear coordinates of nearby edges, the effect of a Whitehead move is illustrated in the 
figure, viz., 



{A, B, C, D, Z) {A + (f>(Z), B - <f>(-Z), C + <f>(Z), D 



-Z),-Z) (21 



In the various cases where the edges are not distinct and identifying an edge with its shear 
coordinate in the obvious notation we have: if A = C , then A' = A + 2(f>(Z); if B = D, 
then B' = B - 2<j>(-Z); if A = B (or C = D), then A' = A + Z (or C = C + Z); if 
A = D (orB = C), then A' = A + Z (or B' = B + Z). 

Sketch of proof Assume that e is the diagonal of a quadrilateral with consecutive sides 
a, b, c, d, where e separates a, b from c, d. Identifying an edge with its lambda length, the 
shear coordinate is given by Z = log ^, i.e., ^ is the required cross-ratio [7j. Further- 
more, if / is the lambda length of the other diagonal, then the lambda lengths satisfy 
Ptolemy's relation ef = ac + bd [7j, and the transformation laws for shear coordinates 
in the proposition are readily derived from this either in the surface for (J2.8)) or in the 
universal cover of the surface in the various cases. q.e.d. 

Insofar as hyperbolic lengths of geodesies are well-defined invariants of homotopy 
classes of curves in F, these lengths must be invariant under Whitehead moves, so we 
have the following 

Lemma 2.4 Transformation ()2.8j) preserves the traces of products over paths 1)2.6)1 . 



2.1.2 Weil— Petersson form 

Th(F) supports its canonical Weil-Petersson Poisson structure, which has a very simple 
form in shear coordinates. 



Theorem 2.5 [22 In the coordinates (Z a 
bracket B WP is given by 



on any fixed spine, the Weil-Petersson 
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dZ Vi dZ Vi+l 



(2.9) 



where the sum is taken over all vertices v and Vi, i = 1, 2, 3 mod 3 ; are the labels of the 
cyclically ordered half- edges incident on this vertex. 

The proof [22] relies on the independence of this form under Whitehead moves as in 
|10j . Indeed, the equivalent expression for the Weil-Petersson Kahler two-form in the 
punctured case was first given in lambda length coordinates in starting from Wolpert's 
formula |HZ| , and the formula in shear coordinates follows from direct calculation using the 
expression Z = log — . (There is more to this geometrically, however, and one must show 
[T§] that the same expression is the Weil-Petersson form for surfaces with boundary.) 

The set of Casimir functions is described by the following proposition, whose proof is 
given in Appendix B. 

Proposition 2.6 The center of the Poisson algebra ()2.9j) is generated by elements of the 
form J2Z a , where the sum is over all edges of T in a boundary component of F{T) and 
the sum is taken with multiplicity. 

2.1.3 Mapping class group description using graphs. 

Recall that the mapping class group MC(F) of an open surface F is the group of ho- 
motopy classes of orientation-preserving homeomorphisms of F. No special constraints 
are imposed by the circle boundary components, i.e., a homeomorphism must fix each 
boundary component only setwise, and the homotopies must likewise fix each boundary 
component only setwise. Thus, if F has b boundary component circles, p punctures, and 
genus g, then MC(F) MC(F g +p ), so we generally write MC g = MC(F) for any surface 
of genus g with s boundary components. In this section, we establish the combinatorial 
presentation of MC g associated with the cell decomposition of decorated Teichmuller 
space. 

Recall that a cell in the decomposition of T g is described by the homotopy classes of 
an embedding i : V — > F of a fatgraph T as a spine of F. MC s g acts on the set of homotopy 
classes of such embeddings by post-composition, and the cell decomposition of decorated 
Teichmuller space T g descends to an orbifold cell decomposition of M. s g = T g /MC(F). 

The modular groupoid MG s g = MG(F) is the fundamental path groupoid of M. g , and 
MC g arises as the subgroup of paths based at any point. Specifically, consider the dual 

graph Q s = Q{F) of the codimension-two skeleton of this decomposition of M. s g (where 
there is one vertex for each top-dimensional cell, edges correspond to Whitehead moves, 
and two-dimensional cells correspond to pairs of homotopic paths in the one-skeleton 
which are homotopic to real endpoints in M s g .) The fundamental path groupoid of Q g is 
the modular groupoid, and in particular, MC g is the stabilizer in MG g of any vertex of 

Go- 

We may think of a Whitehead move along edge a of fatgraph V producing another 
fatgraph T a as an ordered pair (r, T a ), i.e., an oriented edge of Q s . Letting [r 1; T 2 ] denote 
the MC*-orbit of a pair (r!,r 2 ) by the diagonal action of the mapping class group, the 
natural composition descends to a well-defined product 



[r!,r 3 ] = [r ls r 2 ][r 2 ,r 3 ]. 



Theorem 2.7 The modular groupoid MG s g is generated by Whitehead moves and 

relabelings by fatgraph symmetries. A complete list of relations in MG s g is given by rela- 
belings under fatgraph symmetries together with the two following relations. 

Commutativity If a and (3 are two edges with no common endpoints, then 

[r Q( g,r a ][r a ,r] = [r^r^r^r]. 

Pentagon If a and (3 share exactly one common endpoint, then (see Figure 2 drawn 
for the dual graph) 

[r, r a ] [T a , Tp a ] \Tj3a, T a p] [T a p, Tp] [Tp, T] = 1. 

Furthermore, the expression in Theorem \2.5\ for the Weil-Petersson form is invariant 
under Whitehead moves. 




Figure 2-pentagon identity 

Proof of Theorem The first parts are immediate consequences of the cell decomposition. 
Specifically, by connectivity of T g s , any two points can be joined by a smooth path, which 
we may put into general position with respect to the codimension-one^faces; this proves the 
first part. For the second part, a homotopy between edge-paths in A4 s g can likewise be put 
into general position with respect to the co dimension-two faces; there are two possibilities 
for a pair of edges depending upon whether their vertices are disjoint or not, corresponding 
respectively to the commutativity and pentagon relations, proving the second part. The 
invariance of the expression for the Weil-Petersson form under Whitehead moves is a 
direct calculation in lambda lengths ^T] or shear coordinates using (j2.8|) . q.e.d. 



2.2 Poisson algebras of geodesic functions 

The algebra generated (with multiplication and with the Weil-Petersson Poisson bracket) 
by the functions {G 7 } ()2.6|) was first studied by W. Goldman [2"3~j . 

2.2.1 Multicurves 

In the sequel, disjointly embedded families of geodesies will play a special role as they 
constitute a basis for the algebra of observables in both the classical case considered here 
and the quantum case discussed in Section 3. The homotopy class of such a family is 
called a multiple curve. A multicurve is multiset based on the set of curves in a multiple 
curve. 

Definition 2.1 Consider the homotopy class of a finite collection C = {71, . . . , 7„} of 
disjointly embedded (unoriented) simple closed curves 7$ in a topological surface F, where 
C need not be a mutiple curve. A generalized multicurve (GMC) C in F is a multiset 
based on C; one thus imagines Sj > 1 parallel copies of components of C, or in other 
words, positive integral weights s, on each component of C, where Sj is the multiplicity 
of 7j in C. Further, given a hyperbolic structure on F, we associate to C the product 
Gq = Gf* . . . Gf? of geodesic operators ()2.6|) of all geodesies constituting a GMC; these 
operators Poisson commute in the classical case since the components of C are disjoint. 
In particular, a GMC containing a contractible component (of length zero) is twice the 
GMC with this curve removed. 

An edge-path on a spine r C For its corresponding geodesic 7 in F is said to be 
graph simple with respect to Y if it does not pass more than once through any edge of Y. 
Obviously, the set of graph simple geodesies depends upon Y and is not invariant under 
Whitehead moves. Nevertheless, this notion will be useful in what follows. 

2.2.2 Classical skein relation 

The trace relation tr (AB) + tr (AB^ 1 ) — tr A- tvB — for arbitrary 2x2 matrices A and B 
with unit determinant allows one to "disentangle" any product of geodesic functions, i.e., 
express it uniquely as a finite linear combination of GMCs. Introducing the additional 
factor j^-G to be the total number of components in a GMC, we can then uniformly present 
the classical skein relation as 




(2.10) 



2.2.3 Poisson brackets for geodesic functions. 

Turning attention now to the Poisson structure, two geodesic functions Poisson commute 
if the underlying geodesies are disjointly embedded. By the Leibnitz rule for the Pois- 
son bracket, it suffices to consider only "simple" intersections of pairs of geodesies with 



respective geodesic functions G\ and G 2 of the form 

d = ti l ...X l c R l X l z L l X\... ) (2.11) 
G 2 = tr 2 . . . X 2 B L 2 X 2 Z R 2 X 2 D . . . , (2.12) 

where the superscripts 1 and 2 pertain to operators and traces in two different matrix 
spaces. 

The bracket between Xq and X 2 B possesses a simple r-matrix structure 

{Xh,X 2 B } = \{-lTK4i®el)Xh®Xl, (2.13) 

where the "elementary" matrix has entry unity in its ith row and jth column and zero 
otherwise. Direct calculations then give 

{G 1 ,G? 2 } = i(G H -G I ), (2-14) 

where G Y corresponds to the geodesic that is obtained by erasing the edge Z and joining 
together the edges 11 A n and "£)" as well as a B v and "C" in a natural way as illustrated 
in the middle diagram in (|2.10|) : G H corresponds to the geodesic that passes over the 

edge Z twice, so it has the form tr . . . XqRzRd XbLzLa ... as illustrated in the 

rightmost diagram in (j2.1()jl . These relations were first obtained in in the continuous 
parametrization (the classical Turaev-Viro algebra). 

Torus Example For the torus, Th(F^) has three generators X, Y, Z, where 

{X,Y} = {Y,Z} = {Z,X} = 2, 

corresponding to the combinatorially unique cubic spine and the Casimir element is X + 
Y + Z. The geodesic functions for the three graph simple geodesies are 



G x = tr LX Y RX Z = 


e ' 


-Y/2- 


-Z/2 


+ e 


-y/2-t 


-Z/2 _j_ e Y/2+Z/2 


G Y = tr RX X LX Z = 


e ' 


-Z/2- 


-X/2 


+ e 


-Z/2- 


,X/2 + e Z/2+X/2\ (215) 


G z = tr RX Y LX X = 


e ' 


-X/2- 


-Y/2 


+ e 


-X/2- 


fy/2 + e X/2+Y/2_ 


Introducing the geodesic function 














Gz = tr L Y R Z R X L Z = e ^ x/2 - Y/2 ' 


-z 


+ e 


X/2- 


-y/2/ 


e~ z 


+ e Z + 2)+ e X/2+Y/2 + Z^ 



obtained from Gz by a Whitehead move, we find that {Gx,Gy} = Gz/2 — Gz/2, and 
because relation ()2.1U|) implies that GxGy = Gz + Gz, we have 

{Gx, Gy} = -GxGy — Gz, (2.16) 

plus the cyclic permutations in X, Y, Z, i.e., the classical Poisson algebra closes in the 
algebraic span of the geodesic functions {Gx, Gy, Gz}- 



2.2.4 Poisson geodesic algebras for higher genera 



In order to generalize the torus example, we must find a graph on which graph simple 
geodesies constitute a convenient algebraic basis. Such a graph is illustrated in Fig- 
ure 3, where m edges pairwise connect two horizontal line segments. Graph simple closed 
geodesies in this picture are those and only those that pass through exactly two differ- 
ent "vertical" edges, and they are therefore enumerated by ordered pairs of edges; we 
denote the corresponding geodesic functions Gij where i < j. The Poisson algebra for the 
functions Gij is described by 



{Gij, Gkl} 



0, j < k, 

0, k < i, j < I, 

GikQji -QkjQu, i < k < j < l, 

— Ga, j = k, 
Gil - \GijGu, i = k, j <l 
Gik — \GijGkj, j = I, i < k. 



(2.17) 



The graph in Figure 3 has genus y — 1 and two faces (holes) if m is even and genus 
(m— 1)/2 and one face (hole) if m is odd. Such geodesic bases for m even were considered 
in [21]. The Poisson algebras of geodesies obtained there coincide exactly with (|2.17|) . 
These are the so-called so q {m) algebras whose representations were constructed in [2*5] . 

In the mathematical literature, this algebra has also appeared as the Poisson algebra 
of the monodromy data (Stokes matrices) of certain matrix differential equations 
and on the symplectic groupoid of upper-triangular matrices G |2Zj • These matrices have 
entries given by unity on the main diagonal (i.e., we set Gu = 1) and the entries Gij above 
it. For m x m- matrices, there are [y] central elements of this algebra generated by the 
polynomial invariants fg(X) = det(£? + \G T ) = Yl fi{G)X l - The total Poisson dimension 
d of algebra (J2~T7|) is m( m ~ 1] - [f ] , and for m = 3, 4, 5, 6, ... we have d = 2, 4, 8, 12, . . . . 
The dimensions of the corresponding Teichmuller spaces are D = 2, 4, 8, 10, . . . , so we see 
that the Teichmuller spaces are embedded as the Poisson leaves in the algebra ()2.17|) . 




J ■■■ i 

Figure 3-the special spine for higher genera 
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3 Quantization 



A quantization of a Poisson manifold, which is equivariant under the action of a discrete 
group T>, is a family of *-algebras A h depending on a positive real parameter h with D 
acting by outer automorphisms and having the following properties: 

1. (Flatness.) All algebras are isomorphic (noncanonically) as linear spaces. 

2. (Correspondence.) For h = 0, the algebra is isomorphic as a "D-module to the 
*-algebra of complex-valued functions on the Poisson manifold. 

3. (Classical Limit.) The Poisson bracket on A given by {01,02} = lim^o 
coincides with the Poisson bracket given by the Poisson structure of the manifold. 

3.1 Quantizing Teichmiiller spaces 

Here we construct a quantization T h (F) of the Teichmiiller space Tjg{F) that is equivariant 
with respect to the action of the mapping class group V = MC (F) . 

Fix a cubic fatgraph T as spine of F, and let T h = T h {T) be the algebra generated by 
Z^, one generator for each unoriented edge a of T, with relations 

[Z*,Z$\=2Trih{z a ,Zf,} (3.1) 

(cf. (12.9)1 ) and the *-structure 

& = Zt (3.2) 

where z a and {•, •} denotes the respective coordinate functions and the Poisson bracket 
on the classical Teichmiiller space. Because of ()2.9|h the righthand side of (|3.1|) is a 
constant taking only five values 0, ±2nih~, and ±4nih depending upon the coincidences 
of endpoints of edges labelled a and (5. 

Lemma 3.1 The center Z h of the algebra T h is generated by the sums ^2 ae j Z^ over 
all edges a & I surrounding a given boundary component, and the Poisson structure is 
non- degenerate on the quotient T h / Z h . 

A standard Darboux-type theorem for non-degenerate Poisson structures then gives 
the following result. 

Corollary 3.2 There is a basis for T h /Z h given by operators Pi, qi, fori = 1, . . . , 6g — 6 + 
2s satisfying the standard commutation relations [pi, qj] = 2^ih5ij, i,j = 1, . . . , 3g — 3 + s. 

Not only is the proof of Lemma 13.11 given in Appendix B, but also an algorithm for 
diagonalizing this Poisson structure is described there. 

Now, define the Hilbert space Tt to be the set of all L? functions in the g-variables 
and let each g-variable act by multiplication and each corresponding p variable act by 
differentiation, pi = 2irih~ For different choices of diagonalization of non- degenerate 
Poisson structures, these Hilbert spaces are canonically isomorphic. 



Torus Example In the case of the bordered torus, we have three generators, X h , Y h , and 
Z h , the commutation relations (jSHJ) have the form [X h } Y h ] = [Y h , Z h ] = [Z h , X h ] = inih, 
and the single central element is X h + Y h + Z h . In the Darboux-type representation, we 
can identify, e.g., (X h + Y h )/2 with q and (-X h + Y h )/2 with (-2nih)d/dq. 

On the level of the modular groupoid as a category where all morphisms are invertible 
and any two objects are related by a morphism, we have constructed one *-algebra per 
object. In order to describe the P-equivariance we must associate a homomorphism 
of the corresponding *-algebras to any morphism in the modular groupoid. For this, 
we associate a morphism of algebras to any Whitehead move and must verify that the 
relations in Theorem 12. 71 are satisfied. 

We now define the quantum Whitehead move or flip along an edge of T by Eq. (J2.8)) 
using the (quantum) function 

4>{z) = <f{z) = -— / — ^—j— dp, (3.3) 
2 J n smh(7rp) smh(ir hp) 

where the contour Q goes along the real axis bypassing the origin from above. The 
function (j3.3J) is Faddeev's generalization JH] of the quantum dilogarithm. 

Proposition 3.3 For each unbounded self-adjoint operator Z on 7i, <p(Z) is a well- 
defined unbounded self-adjoint operator on Ti. 



1+ e~ z ' 
2vri 



Proof The function <p h (Z) satisfies the relations (see j3j) 

<p\z) - <p\-z) = z, 

9irih 

(j) n {Z + iTTh) - <f>\Z - iTTh) 

and is meromorphic in the complex plane with the poles at the points {ni(m+nh), m,n G 
Z + } and {— 7u(m + nh), m,n G 

The function <fi h {Z) is therefore holomorphic in the strip |ImZ| < 7rmin(l,Re^) — e 
for any e > 0, so we need only its asymptotic behavior as Z G ffi. and \Z\ —>■ oo, for which 
we have (see, e.g., |28j) 

<j>\Z) ={Z+\Z\)/2 + 0{l/\Z\). (3.4) 

\Z\^oo 

Therefore, the function (j) h (Z) increases as \Z\ at infinity and represents an operator in Ti. 
by the functional calculus |HH]- q.e.d.. 

Theorem 3.4 The family of algebras T h = T h (T) is a quantization of Th(F) for any 
cubic fatgraph spine T of F, that is: 

1. In the limit h i— ► 0, morphism ()2.8|) using \3. 5j) coincides with the classical 
morphism ()2.8|) where <f)(Z) = log(l + e z ); 



2. Morphism ()2.8|) using (|3.3|) is indeed a morphism of *- algebras; 

3. v4 /Zip satisfies W% = I, cf. (|2.8|) . and flips satisfy the commutativity relation; 

4. F/zps satisfy the pentagon relation. 1 
Furthermore, 

5. TTie morphisms T h (T) —>■ T l / n (Y) given 6?/ \— > Z& commute with morphisms 

(ED. 

Sketch of the proof Property 1 follows since lim^ 4> h {z) = log(e 2 + 1), and Property 3 
is obvious. 

In order to prove Property 2, we must first verify that [A + <ft h (Z), B — <j) h (—Z)] = 
and [A + 4> h (Z), D — (p h (— Z)\ = —2irih~ (since the other relations are obviously satisfied), 
which follows from the identity <fi h (z) — <fi h (—z) = z. 

For Property 5, we must verify that the morphism T h (T) — > T^^r) commutes with a 
flip, that is, (A + (j) h (Z))/h = A/h + (f) h (Z/h), (B - <j) l l h {-Z))/h = B/H- <jfi(-Z/h), etc., 
which follows from (j) h (z)/h = ^^(z/h). 

Turning finally to the most nontrivial Property 4, we may reformulate it as follows. 
There are seven generators involved in the sequence of flips depicted in Figure 2 for the 
dual cell decomposition, which are denoted A, B, C, D, E, X, Y as in the figure. As a 
result of a flip, the piece of graph shown in Figure 1 just gets cyclically rotated. Denote 
by Aj, Bi,Ci,Di,Ei,Xi, and Yi the algebra elements associated to the edges of this piece 
of graph after % flips are performed. From (|2.8|) . (|3.3|) . these elements evolve as follows: 



/ Xi+i 










\ 


Yi+1 






-x, t 






A i+1 






A 






Bi + i 






Ei 






Ci+l 






A + 






A+i 






Bi — 


i> h {-Xi) 




\ E i+ i 


) 




\ Ci + i 




J 



We must prove that this operator is periodic with period five. 

Assume for a moment that this five-periodicity of X$ has been established. Then 
five-periodicity of Yi and other variables follow from simple calculations. 

It suffices to prove this five-periodicity for the Xi since the five-periodicity of the other 
operators such as Y$ then follows from elementary calculations. As to the five-periodicity 
of Xi, let us "take logarithms" and introduce four new algebra elements 

Ui = e x >; V, = e Y >- Ui = e x ^ h ; V, = e Y >/\ 

which satisfy the following commutation relations 

UtVi = q~ 2 ViU i} UiVi = q- 2 ViUu 
UiVi = VjUj, Villi = UiVt, (3.6) 

1 This result was independently obtained by R. M. Kashaev 



where 



q = e ' 7vih , q = e ~ m/h . 



Under the flip, these variables are transformed in an especially simple way, 

Ui-i = Vr 1 (3.7) 

V i - 1 = U i (l + qV i ) (3.8) 

Ur-i = Vr 1 (3.9) 

K_ 1 = y i (l + gK f ). (3.10) 

As the first step of the proof, we consider the inverse transformation laws for Xi and Yf. 

x f _i = -y i; K t _ 1 = x i + /l (^)- 

Equations ()3.7|) and (|3.9|) are obvious. Using the standard formula 

1 (B 
g [A,B] JB 

we obtain 



A+F(B) _ f* +[A > B] F( z )dz A 



y._ } = e Yi-i = e x i +^(Y i ) = e Xi e g +2 " in <)fi(z)dz 

( Ilk f e -ipz l e -2pirh _ ^\ \ / w fl r e -ip(z-nih) 

— XJ ^ I / dp J — XJ % Gxp I / dp 

V 2 Jn (-ip) sinh(7rp) smh.(nhp) J 1 \i J n p smh(irp) 

= U l {l + q- 2 V l ). 

The proof of ()3.10|) is analogous. 

Now, in order to finally prove that Xi is five-periodic, it suffices to verify that both 
Ui and Ui are five-periodic. Indeed, if only the operator Ui is five-periodic, it does not 
suffice because the logarithm of an operator is ambiguously defined. However, if we have 
two families of operators U and U , which depend continuously on H, then, assuming the 
existence of an operator X (depending continuously on h) such that U = e x and U = e x ^ h , 
then this operator is evidently unique. (It can be found as \im^ m+n /f i ^Q(U m U n )/(m+n/h) 
for any irrational value of h.) The five-periodicity of sequence (|3.8|) (and ()3.10|) ) is a direct 
calculation using (|3.6j) . 

The only subtlety remaining is the possibility that some of the edges A, B, C, D, E 
coincide. (Note, however, that X and Y must have exactly one common vertex.) If, say, 
edges A and C coincide, then the value of the commutator [A, X] is doubled by definition, 
and we can then fictitiously split the edge A = C into two half-edges with the matrices 
Xa/2 assigned to each half-edge using the formula 

Xa = X A / 2 ( ) Xa/2- (3-11' 




The quantities A/2 = C/2 have the same commutation relations with the rest of variables 
as well as the same transformation laws as the quantities A and C before, so the earlier 



formulas remain valid if we simply replace there A' by A'/ 2 and Aby A/2. The net effect 
is that commutators with A = C are doubled. 

If edges A and B coincide, then we must correct formulas (|3.4|) using a splitting as 
above but demanding A'/2 = A/2 + X/2. Obviously, [A, X] = in this case (in which 
formulas for the quantum ordering below must be also corrected, see Section 3.5). 

The formulas in Proposition 12.31 can thus be realized for exponentiated quantities, 

although in the current quantum case there will be corrections. For instance, formulas 

(|3.8|) and ()3.9|) will be different; indeed, letting A = C = X, we calculate that e x ^ 2 i— > 
e x/2+<f> h (z) = . e x/2(i _|_ e Z) .= e x/2 _|_ e x/2+z = e x/2^ + e -™h e z^ gince the commu t a tor 

of X and Z is doubled, where the normal ordering (the Weyl ordering) : • : is explained 
in the next section. The transformation for e x itself becomes more complicated when 
A = C, namely, e x = (e x/2 ) 2 i-> e x (l + q 3 e z )(l + qe z ). q.e.d. 

Corollary 3.5 1. Let K be an operator acting in the Hilbert space L 2 (IR) and having the 
integral kernel 

K(x,z) = F h (z)e~^, (3.12) 

where 



F\z) = exp -- / — - ; — - dp (3.13) 
\ 4 J n psmh{irp) smh{7r hp) J 

Then the operator K is unitary up to a multiplicative constant and satisfies the identity 

K 5 = const. (3.14) 

2. Let h = m/n be a rational number and assume that both m and n are odd. Introduce a 
linear operator L(u) acting in the space C n and depending on one positive real parameter 
u through its matrix 

L(u)) = F\ 3 ,u)q-^ } (3.15) 

where 

F\j, u) = (1 + u) j ' n JJ(1 + q- ik+2 u^ n )-\ 

k=0 

Then the following identity holds: 

L(u)L{v + uv)L(v + vu- 1 + u'^Liu^v' 1 + u -1 )^ -1 ) = I. (3.16) 



Using this construction, Kashaev [28] constructed the set of eigenfunctions of the 
quantum Dehn twist transformation. Namely, Kashaev's dilogarithm function e^z) is 
1/F h (2z) from ()3.13|) . We return to this discussion when considering sets of quantum 
Dehn twists for the torus in Section 5. 

3.2 Geodesic length operators 

We next embed the algebra of geodesies (j2.6J) into a suitable completion of the constructed 
algebra T h . For any 7, the geodesic function G 7 (12.6)1 can be expressed in terms of shear 



coordinates on T H : 



tTP Zl -z n = Xl ex P \ \ E m j(li a ) z * \ , (3-17) 



where ^ (7, ct) are integers and J is a finite set of indices. In order to find the quantum 
analogues of these functions, we denote by T h a completion of the algebra T h containing 
e xZa for any real x. 

For any closed path 7 on F, define the quantum geodesic operator GO: G T h to be 



G*=*trP Zl ... Zn 



7 



= E ex P \l E K(7,«)^+ 2^(7, a)) 1 , (3.18) 



where the quantum ordering * ■ * implies that we vary the classical expression ()3.17|) 

by introducing additional integer coefficients 0^(7, a), which must be determined from 

1 

the conditions below. Notice that the operators {G7} themselves can be considered as 
belonging to the algebra T h insofar as 



G* = ^exp{— ^ (m j ( 1 ,a)Z2 + 2nic^,a)) }. (3.19) 

In what follows for the notational simplicity, we shall sometimes omit the superscript 
h from G h and write it merely G assuming that G is either an operator or a classical 
geodesic function depending on the context. In this paper, we shall concentrate on the 
case of quantum functions of the G h -type; the consideration of the sector G x ' h is analogous 
and does not lead to new effects (at least at the present stage of understanding), so we 
omit it. We also call a quantum geodesic function merely a quantum geodesic (since 
implicitly quantum objects admit only functional, not geometrical, descriptions). 

We wish to associate an operatorial quantum multicurve QMC to a multiset C of 
quantum geodesies corresponding to disjointly embedded families of nonnegative integrally 
weighted geodesies. One ansatz will be that operators corresponding to disjoint underlying 
geodesies must commute; this implies that the ordering in which the quantum geodesies 
enter the product QMC is immaterial, where the product is defined as for GMCs. We 
next formulate the defining properties of quantum geodesies. 

1. If closed paths 7 and 7' do not intersect, then the operators G^ and Gz, commute. 

2. Naturality. The mapping class group MC(F) (|2.8j) acts naturally, i.e., for any 
{G^}, 5 G MC(F) and closed path 7 in a spine T of F, we have 8(G*) = Gf r 

3. Geodesic algebra. The product of two quantum geodesies is a linear combination 
of QMC's governed by the (quantum) skein relation [29J. 



4. Orientation invariance. Quantum traces of direct and inverse geodesic operators 
coincide. 



5. Exponents of geodesies. A quantum geodesic G„ 7 corresponding to the n-fold 
concatenation of 7 is expressed via G 7 exactly as in the classical case, namely, 

G nj = 2T n (G 7 /2), (3.20) 

where T n (x) are Chebyshev's polynomials. 

1 

6. Duality. For any 7 and 7', the operators and G*, commute. 

We shall let the standard normal ordering symbol :e ai e a2 ■ ■ ■ e a ": denote the Weyl 
ordering e aiH ha ", i.e., 

:e ai e a2 ■ ■ ■ e a ": = 1 + (01 H h a n ) + j {a x + • •• + a n )(a 1 H + a n ) H (3.21) 

for any set of exponents with a, 7^ —aj for i ^ j, In particular, the Weyl ordering implies 
total symmetrization in the susbscripts. 

Proposition 3.6 For any graph simple geodesic G with respect to any spine T, the coef- 
ficients c*(7, a) in (j3.18j) are identically zero. 

Proof Consider term-by-term the trace of the matrix product for the quantum graph 
simple geodesic and expand it in Laurent monomials in e Zi ^ 2 . It is easy to see that 
each term e Zi l 2 comes either in power +1, or —1 in the corresponding monomial and 
there are no equivalent monomials in the sum. This means that in order to have a 
Hermitian operator, we must apply the Weyl ordering with no additional g-factors (by 
the correspondence principle, each such factor must be again a Laurent monomial in q 
standing by the corresponding term, which breaks the self-adjointness unless all such 
monomials are unity). Since quantum Whitehead moves must preserve the property of 
being Hermitian, if a graph-simple geodesic transforms to another graph-simple geodesic, 
then a Weyl-ordered expression transforms to a Weyl-ordered expression, and only these 
expressions are self-adjoint. q.e.d. 

Torus Example For the torus with one hole, there are three graph simple quantum 
geodesies for any spine, which are exactly ()2.15|) in the Weyl-ordered form. The quantum 
geodesies Gz obtained from Gz by the flip transformation is 

G z = e - x / 2 ~ y / 2 - z +e x / 2 - y / 2 - z +e x / 2 -^ (3.22) 

The product of two graph simple quantum geodesies is 

G X G Y = e i7rh/2 G z + e~ inh/2 G z . (3.23) 

Denoting q = e ~ mf \ [A, B] q = q x l 2 AB — q~ l l 2 BA, and £ = q — q -1 , we obtain from 

[Gx, Gy] q = £,Gz, [Gy, Gz] q = £,Gx, [Gz,G x ] q = ^Gy- (3.24) 

This algebra is exactly the so 5 (3) quantum algebra studied in [23 • There is a unique 
central element, the quantum Markov relation 

M = GxGyGz - q 1/2 (G 2 x + q- 2 G 2 Y + G 2 Z ). (3.25) 



3.3 Algebra of quantum geodesies 



Let G± and G2 correspond to the respective graph simple geodesies with respect to the 
same spine having one nontrivial intersection. For G\ and G2, formula ()3.18|) implies, by 
virtue of Proposition 13.61 the mere Weyl ordering. 
After some algebra, we obtain (cf. (|2.14j) ) 

G X G 2 = e~ i7rh/2 G z + e i7rh/2 G z , (3.26) 

where Gz coincides with the Weyl-ordered G I in the classical case (cf. (12.14(0 while Gz 
contains the quantum correction term 

Gz = -txHx 2 ...{e\ j ®e%)[X 1 z ®Xl]...l 

= :tr 1 tr 2 ...(eJ i ®ep[X^(g)X| + 2(l-cos7rn)e} 1 ®ey ... : . 

Here eh (8> is the standard r-matrix that permutes the spaces "1" and "2," and as a 
result, the "skein" relation of form (j2.14|) appears. Locally, this relation has exactly the 
form proposed by Turaev [2H|, i.e., for two graph simple geodesies intersecting at a single 
point, we have the defining relation 

Gz Gz 



' 1 (3.27) 

(The order of crossing lines corresponding to G\ and Gi depends on which quantum 
geodesic occupies the first place in the product; the rest of the graph remains unchanged 
for all items in ()3.27j) ). Note, however, that if the quantum geodesies G\ and G2 cor- 
respond to graph simple geodesies, we may turn the geodesic Gz again into the graph 
simple geodesic G' z by performing the quantum flip with respect to the edge Z. 

If we now compare two unambiguously determined expressions: G' z , which must be 
Weyl ordered, and Gz obtained from the geodesic algebra, we find that Gz = G' z . This 
enables us to formulate the main assertion. 

Lemma 3.7 j3] There exists a unique quantum ordering * . . . * (I3.18J1 . which is generated 
by the quantum geodesic algebra \3.21 ) and consistent with the quantum mapping class 





groupoid transformations (|2.8|) . i.e., so that the quantum geodesic algebra is invariant 
under the action of the quantum mapping class groupoid. 

We can now relax the constraints of graph simplicity of curves: as the quantum 
geodesic algebra is quantum mapping class group invariant, having two arbitrary embed- 
ded geodesies with a single intersection, we can transform them using quantum morphisms 
to a canonical form of graph simple geodesies and employ the Weyl order. Relation ()3.27|) 
remains valid in both cases. 

Let us now address the problem of multiple intersections. Here, we have the following 
lemma. 



Lemma 3.8 [4 If more than one intersection of two QMCs occurs, then the quantum 
skein relations (|3.27|) must be applied simultaneously at all intersection points. 

This lemma implies the standard Reidemeister moves for curves on a graph, where the 
empty loop gives rise to a factor — e ~ lnh — e mh \ that is, for geodesies intersecting gener- 
ically, apply ()3.27|) simultaneously at all intersection points to obtain the Reidemeister 
moves. 

Remark 3.1 The quantum algebra Ail was studied in |4j, where the exact correspon- 
dence with Kauffman bracket skein quantization of the corresponding Poisson algebra of 
geodesies (see (201) was observed. In j2Z], this algebra arises as the quantum deformation 
algebra of the classical groupoid Poisson relations for the group SX(4). This algebra is 
however only one among many quantum Nelson-Regge algebras corresponding to Rie- 
mann surfaces of higher genera which are described in the next section. 



3.4 Quantizing the Nelson— Regge algebras 



The algebra (|2.17|) was quantized by the deformation quantization method in (2UEH]. We 
now explicitly implement the quantization conditions (|3.1|) . It is convenient to represent 
the (classical or quantum) elements Qij as chords connecting the points of the cyclically 
ordered set of indices There is then a trichotomy: if two chords do not intersect, 
then the corresponding geodesies do not intersect either, and the quantum geodesies 
commute (Figure 4a); if two chords have a common endpoint, then the corresponding 
geodesies intersect at one point, and the three quantum geodesies Qfp Qj k , Q\ { (as depicted 
in Figure 4b) constitute the quantum subalgebra so q (3); if two chords intersect at an 
interior point (as depicted in Figure 4c), then the corresponding geodesies intersect in 
two points, and the corresponding quantum geodesies and Q kl , i < k < j < I, satisfy 
the commutation relation 

$ki\ = ZiGikGji ~ & h jk) ( 3 - 28 ) 
with the usual commutator (not the g-commutator) and where again £ = q — q^ 1 . 






(3.29) 



Figure 4-picture of Nelson— Regge quantum relations 



3.5 Improving the quantum ordering 

We now extend the construction of QMCs by taking the products of operatorial matrices 
Xz, L, and R along oriented geodesies as before, but we do not apply the trace operation 
as in the definition of the geodesic operators. 



For the three cases of oriented curves depicted in Figure 5 below, we may apply the 
indicated quantum Whitehead move and calculate as follows: 



For curve 1, we obtain 

e -B'/2 g -Z/2 g A'/2 _|_ g -B'/2 g Z/2 g A'/2 _ g -B'/2 g Z/2 g -A'/2 
& B> /2 & -Z/2 & A' /2 

e A/2-B/2 _ g -A/2-B/2- J [A,B] 

e A/2+B/2-i[A,B] q 



for 5' ^ A' 



/ -B/2 A/2 -B/2 -A/2 \ 

e-^i 66 e e for B' ^ A' 

\ e B ' 2 e A ' 2 J 

r 3.30 
X B LX A for A = B. 



For curve 2, we obtain 
X(ji RX_zLXa' = 



e -C/2 g Z/2 g A'/2 _ g -C/2 g Z/2 g -A'/2 

e C"/2 g Z/2 g A'/2 _|_ g C"/2 g -Z/2 g A'/2 _ g C"/2 g Z/2 g -A'/2 

g A/2-C/2+Z/2+i[A,B] _ g -A/2-C/2+Z/2 _ g -A/2-C/2-Z/2 
e A/2+C/2+Z/2 _ g - A/2+C/2+Z/2- i [A,B] 

e -C/2 g Z/2 g A/2 _ g -C/2 g Z/2 g -A/2 _ g -C/2 g -Z/2 g -A/2 
e C/2 g Z/2 g A/2 _ g C/2 g Z/2 g -A/2 

= X C LX Z RX A , (3.31) 

and if edges A and C coincide, we merely use the same splitting as in (|3.11j) : no additional 

factors arise. 

For curve 3, we obtain 



Xqi RX—Z RXjji 



_g-C'/2gZ/2g-D'/2 
e C/2 g -Z/2 g C/2 _ g C/2 g Z/2 g -C/2 _ g -C/2 g -Z/2 g -C/2 



_ e -C/2-D/2+i[A„B] 
e C/2+D/2+i[A,B] _ e C/2-D/2 



for c ^ jy 



( -p-C/2 -D/2 \ 

el^M , , , 6 , forC'^D' 

e i[ A ' B lX Cj RX D . forC^D, 
X C RX D for C = D. 

In these formulas, we have used the identities 

e A'/2+C/2+Z/2 + e A'/2+C/2-Z/2 = . g A/2+C/2 1 / g Z/2 j g -Z/2> 

1 + e ~ z 



(3.32) 



_ e A/2+C/2+Z/2 



(3.33) 



and 

e -A'/2~C'/2+Z/2 = . e -A/2-C/2/ 1+ e -Z\ e Z/2 . 

_ e -A/2~C/2+Z/2 + e -A/2-C/2-Z/2^ (3M) 

where vertical dots denote the Weyl ordering (j3.21|) as before implying total symmetriza- 
tion with respect to all the variables {A, B, C, D, Z}. 

As we have just observed and amazingly enough, the only thing that changes under a 
quantum Whitehead move is the overall factor standing by the product of matrices, and 
even this factor can easily be taken into account if included with each left-turn matrix 
L is an overall factor e~^ A,B \ and included with each right-turn matrix R is an overall 
factor e a^ A ' B K 

The subtlety of potentially multiple intersections between curves corresponds to the 
possibility that edges of the graph may coincide. As in the previously useful convenient 
fiction of splitting such an edge into two half-edges, the formulas ()3.30|) . ([3.31)1 . and ()3.32|) 
remain valid if we replace there A'/ 2 by A' /A and A/2 by A/4 in case A = C for instance. 

We thus have the quantum analogue of Lemma 12.41 

Lemma 3.9 For any oriented non boundary-parallel geodesic 7, taking the finite (peri- 
odic) sequence of matrices with quantum entries in X Zi as before but making the replace- 
ment 

L = q~ 1/4 L, R = q 1/4 R, (3.35) 
where q = e ~ lnh ! the resulting product of matrices 

P 7 = {X z lX Zn _ x R . . . Xz 3 LX Z2 RX Zl R) (3.36) 

is invariant under quantum Whitehead moves. 

Remark 3.2 In the case of a boundary-parallel curve, where a sub-word of the form 
Xj)iRX_zR transforms to a sub- word of the form XpR, we must set the resulting turn 
matrices to be 

L = q~ 1/2 L, R = q 1/2 R, 
i.e., they have quantum factors doubled in comparison with L and R. 

We now address the question of the proper quantum ordering. Using MCG transfor- 
mations, we can reduce any simple curve either to the form of a graph simple curve with 
exactly one left and one right turn (c.f. formulas (|2.15j) and Figure 5) if this curve is not 

boundary-parallel, or to the form XpR or XpL for a boundary-parallel curve. In both 
cases, each term in the corresponding Laurent polynomial must be self-adjoint, which 
immediately results in the Weyl ordering by Proposition 13. (il 




Figure 5-three cases of flips for geodesies 

This observation does not suffice to derive the proper quantum ordering in our quanti- 
zation of Thurston theory since we cannot consider in a consistent way an infinite product 
of matrices corresponding to infinite leaves of a measured foliation. Our tools for analyzing 
the Thurston theory in this case will devolve to naturality of the MCG action on the QMC 
algebra and an operatorial version of infinite continued fractions in Section 15.3.41 On the 
other hand, the improved quantum ordering is used in our analysis of closed geodesies on 
the torus, i.e., of operatorial finite continued fraction expansions in Section IB. 3. II 



4 Classical Thurston Theory of Surfaces 



Let Fg denote an oriented smooth surface with s > punctures (so F* may be closed 
without boundary in this section), with genus g > 0, and with negative Euler characteristic 
2-2g-s < 0. 

4.1 Measured foliations and Thurston's boundary 

Define a measured foliation on F = Fg to be a one-dimensional topological foliation T of 
F, where in a neighborhood of any p G F, T must restrict (in an appropriate chart) to the 
horizontal foliation as in Figure 6a or it must restrict to a foliation with one n-pronged 
singularity at p, for n > 3, as illustrated in Figure 6b. 




a b 
Figure 6-pictures of foliations 



Furthermore, T comes equipped with a transverse measure p, which assigns to any arc a 
in F that is transverse to T a real number p{a) e M>o, where p is required to satisfy: 

No Holonomy: If a , a± are homotopic through arcs a t transverse to J 7 , for < t < 1, 
keeping the endpoints of a t on the same leaf for all t, then p(a ) = p(ai); 

a-Additivity: If a is the serial concatenation of transverse arcs a±, a%, . . ., then /i(a) = 
Ei>iMai)- 

In other words, in the neighborhood of a non-singular point, there is a local chart 
: U — > M 2 = {(x, y) : x,y & K} so that _1 ({?/ = constant}) are the leaves of the 
foliation J 7 nU. If two charts C/j and C/j intersect, then the transition function 0^ are of 
the form (f>ij(x,y) = (hij(x,y),Cij ± y), where c^- is constant. In these coordinates, the 
transverse measure is \dy\. In case the transition functions can be chosen with constant 
sign <pij(x,y) = (hij(x,y),c + y), i.e., if the foliation is "transversely orientable", then 
(away from the singular points) y is the primitive of a closed one-form on F. 

Another canonical construction of a measured foliation on a Riemann surface is given 
by taking the leaves of the foliation to be the level sets of a harmonic function, where the 
transverse measure is given by integrating the conjugate differential along transverse arcs. 
Still another example is given by the homotopy class of a finite collection of disjointly 
embedded (weighted) curves, as we shall see. 



Figure 7- Whitehead collapses on foliations 



There is an equivalence relation on measured foliations generated by isotopy and 
Whitehead collapse as illustrated in Figure 7, and the set of all equivalence classes (in- 
cluding the empty measured foliation 0) is denoted M.J-' = M.J-'(F), where the class 
of (J-,[j) is denoted [J 7 , fj]. To naturally topologize M.J 7 , we introduce the discrete set 
S = S(F) consisting of all free homotopy classes [c] of simple closed curves c in F which 
are neither null homotopic nor puncture-parallel. We shall also require the set S' = S'(F) 
consisting of homotopy classes of all disjointly embedded families of curves in F, where 
each component of the family lies in S and no two components are homotopic. 

One can show that for any [c] G S, there is a representative Cjr G [c] which 
minimizes the /x-transverse measure in its homotopy class; furthermore, given equivalent 
measured foliations (J r i ,iii), for i = 1,2, we have /^(c^J = /z 2 (cjr 2 ) j^S], and hence there 
is a well-defined mapping 

J : MT -> Rf 
(F,pi) ^ : [c] ^ nipjr) J , 

where the empty measured foliation is identified with = {0} 5 . This mapping J is 
an injection [^nj and induces a topology on M.T (where a neighborhood of is homeo- 
morphic to a cone from over AiJ r /M, >0 , with the natural action by homothety of M>o 
on measures). The function i(r,n) is called the (geometric) intersection function of the 
measured foliation (J 7 ,^). 

Given [c] G S, there is a corresponding measured foliation defined as follows. Choose 
a representative c of [c] and collapse F — c onto a spine (making further choices) to build 
a foliation T c of F, whose leaves either lie in the spine or are homotopic to c, and choose 
a transverse measure \i c on T c that pulls back under the collapsing map the counting 
measure 5 C on c. It is a classical fact due to Whitehead that the resulting Whitehead 
equivalence class [J-" c , /i c ] is well-defined independent of any choices. Taking the projective 
class of this foliation, we may thus regard 

S C VTq. 

Furthermore by construction, i[^ c ,fi c ]([d]), for [d] G S, is just the geometric intersection 
number of [c] and [d], i.e., the total number of intersections of representatives c with d, 



where c and d intersect minimally. More generally, given a family of curves c\ , . . . , c n 
representing a point of S', together with a collection wi,...,w n G M>o of "weights", 
we may again collapse to a spine of F — U{q}™ to produce a foliation T and choose a 
measure p on J 7 that pulls back the weighted sum ^[ Wi5 Ci to get a well-defined Whitehead 
equivalence class [J 7 , p] . As mentioned before, one may thus associate a measured foliation 
to a weighted curve family. 

Of special interest is the subspace M.J-'q = M.J-'q(F) consisting of all measured fo- 
liations p) of compact support, i.e., any leaf of T with a transverse arc a so that 
/i(a) > must be disjoint from a neighborhood of the punctures and no such leaf is 
puncture parallel. There is a natural M>o-action on M.T — {0} and — {0} given 

by scaling the transverse measure, and the corresponding quotients 

VT = VT{F) = {MT{F) - {0})/R >o , 

VF Q = PF (F) = {M!Fq{F) - {0})/M >o 

are the spaces of central interest in the sequel. The projective class of (.F, p) and [J 7 , p], 
respectively, will be denoted (J 7 , p) and [J 7 , p] . 

A point of T = T(Fg) = T g s may be regarded as the class of a hyperbolic metric on F, 

1. e., a complete finite-area Riemannian metric on F of constant Gauss curvature -1. We 
shall also require the "Yamabe space" y = y(F) of all complete finite-area Riemannian 
metrics on F of constant Gauss curvature — x 2 , for some x G M>o- y is canonically 
homeomorphic to T g s x lR>o, where (p, x) corresponds to the class of the metric xp, and 
we let 7r : y — > T denote the projection onto the first factor. Define the map 

/ : y -> Mfo 

where ^ P ([c]) is the p-length of the unique p-geodesic in the homotopy class [c\. Thus, if 
p G y corresponds to (7r(p), x), then I(p) = J(7r(p), x) = x I(tc(p), 1). 

The basic facts [36J are that / : y — > 1R> and J : M-F^F) — > K> are embeddings 
with disjoint images, and we may define a completion y of y in Mx, by setting 

y = i(y) u JiMFv) 

and identifying y with I(y). Passing to quotients under the homothetic actions of M>o 
on R> — {0}, on y — {0}, and on M.J-'o — {0}, we obtain Thurston's compactification 

T B = (T;uPf ) « (3;_ { o})/M >0 
of T « (y - {0})/M >0 by « (A^-^o - {0})/K >o . 

Theorem 4.1 1. 15 , [36J , 32 VJ 7 o{F g ) is naturally a piecewise linear sphere of dimen- 
sion 6g — 7 + 2s which compactifies T g s to produce a closed ball T g . 

2. [H3,[B] The action of the mapping class group MC g on T g extends continuously to 
an action on T s g , where the action on PJF (F|) is the natural one, and there are explicit 
piecewise linear formulas for the action of Dehn twist generators. 



3. Suppose that a sequence of hyperbolic metrics pi on F tends to a point [J-'jp] G 
VTq. In the projectivization of R> — {0}, the projectivized length functions l Pi of p { 
converge to the projectivized intersection function i(F,p) of (J 7 , fi). 

4. [T3], [T7] The function i(T,ix) '■ <S — ► K>o extends continuously to the "geometric inter- 
section pairing" AdJ^o x M.Fo ^>o, that vanishes on the diagonal and which is also 
invariant under MC s g . 

5. [I5].|34|.[19] The Weil-Petersson Kahler two-form on T g continuously extends (in 
the appropriate sense on Yamabe space) to a non-degenerate symplectic form, called 
"Thurston's symplectic form", on A4J r o(F g ) (za y — y), which is invariant under MC s g . 

Though there is all this beautiful natural structure on Thurston's boundary, the quo- 
tient VJ r o{F g )/MCg is maximally non-Hausdorff (i.e., its largest Hausdorff quotient is a 
singleton) as we shall see, so there is no correspondingly nice Thurston compactification 
on the level of Riemann's moduli space. 

Torus Example Recall that for the once-punctured torus F = F^, the Teichmiiller space 
is « {z G C : \z\ < 1}, and the mapping class group is MC s g ~ PSX 2 (Z,). Indeed, the 
right Dehn twists M and L on the meridian and longitude, respectively, generate MC\, 
and a complete list of relations between them is given by i = MLM = LML and l 2 = 1. 

A point of VJ-'o(F) is uniquely determined by its "slope", defined as follows. Fix 
two disjointly embedded ideal arcs x, y asymptotic to the puncture p which decompose F 
into an ideal quadrilateral, where x U {p} is homotopic to the meridian, and y U {p} is 
homotopic to the longitude. Given [J 7 , fj] G A4J-'q(Fi)~{0}, compact support guarantees 
that i[F,n]{x) and i[F,n]{y) are well-defined and finite, and the ratio \9\ = i[r,n]{y) G 
[0, oo] is therefore projectively well-defined; we further imbue 9 with a sign (when it is 
finite and non-zero) in the natural way, where the sign is positive if one (in fact, any) leaf 
of T immediately after meeting x then meets the copy of y in the frontier of the ideal 
quadrilateral which lies to the right (where the orientation of the ideal quadrilateral is 
inherited from that of F±). It is easy to see that the slope 9 is a well-defined and complete 
invariant of [J 7 , p] G VJ r o(F^), where we regard 9 G S 1 in the natural way. Thus, 
VF Q {Fl) w S 1 , and T\ = U TFq(FD is a closed ball (which you should not identify 
with the Poincare disk together with its circle at infinity). The slope 9 = p/q is rational 
if and only if the measured foliation corresponds to the simple closed curve wrapping p 
times around the meridian and q times around the longitude, where p and q are relatively 
prime integers. Geometrically, deforming hyperbolic structure to pinch this curve, it 
is clear from elementary considerations of hyperbolic geometry that the corresponding 
geodesic length functions converge projectively to the geometric intersection number with 
this curve. Furthermore, the geometric intersection number of curves with slopes p/q,r/s 
written in least terms is given by \ps — qr\. 



4.2 Train tracks 



A train track r C F is a graph (where vertices are called "switches" and edges are called 
"branches") together with the following extra structure: 

Smoothness: r is C 1 away from its switches. Furthermore, for each switch v of r, 
there is a tangent line I to r at v in the tangent plane to F at v so that for each 
half-branch whose closure contains v, the one-sided tangent at v lies in £] 

Non-degeneracy: Vertices of r are at least trivalent, and for any switch v of r, there 
is an embedding (0, 1) — > F with /(|) = v which is C l as a map into F; 

Geometry: Suppose that C is a component of F — r, and let -D(C) denote the double 
of C along the C l frontier edges of C so the non-smooth points in the frontier of C 
give rise to punctures of D{C). We demand that the Euler characteristic of D{C) be 
negative. 

The smoothness condition is synonymously called the structure of a "branched one- 
submanifold" and leads to the fundamental notion of a graph smoothly supporting a curve 
or another train track as we shall see. According to the non-degeneracy condition, the 
half-branches incident on a fixed switch decompose canonically into two non-empty sets 
of "incoming" and "outgoing" branches. The geometric condition rules out the following 
complementary regions: smooth disks (i.e., nullgons), monogons, bigons, smooth annuli, 
and once-punctured nullgons, and will be further explained below. 

Let B(t) denote the set of branches of r. A function /i : B(t) — > IR> induces 
ji : {half — branches of r} — > IR> in the natural way (where fi(bi) = fx(b) if bi C b) and 
satisfies the switch conditions provided that for each switch v of r, we have 



Such a function satisfying the switch conditions is called a {transverse) measure on r, 
and r itself is said to be recurrent if it supports a positive measure \i with /i(6) > for 
each branch of r. In the sequel, train tracks will tacitly be assumed to be recurrent. 

Torus Example There is a unique combinatorial type of recurrent trivalent train track 
r in the surface F = F±, and two embeddings of it as spine are illustrated in Figure 8. 
F — t consists of a single once-punctured bigon. There are two branches of r so that /i 
is uniquely determined by its values on these branches, and the weight on the remaining 
branch of r is given by their sum according to the switch condition. Notice that this 
train track is "orientable" in the sense that the graph underlying r admits an orientation 
where incoming points toward outgoing at each vertex. Thus, fixing an orientation on r, 
a measure piir uniquely determines a homology class in Hi(Fl, M). Equivalently, every 
measured foliation of compact support on F± is transversely orientable. 

Construction 4.1 Given a train track r and a positive measure fi on it, we may construct 
a measured foliation of a neighborhood of r in the following way. For each branch b of 



outgoing 
half — branches b 



incoming 
half — branches b 



t, take a rectangle of width /j,(b) and length unity foliated by horizontal leaves. For each 
switch, place the rectangles of the incoming branches next to one another and likewise 
for the outgoing branches, and then finally glue the vertical edges of all the incoming to 
all the outgoing rectangles at each switch in the natural way preserving the transverse 
measure along the widths of the rectangles by the switch conditions. This produces from 
jj, a measured foliation of a tie neighborhood of r, where a vertical leaf in any rectangle 
is called a tie, and the singular ties arise from the vertical sides of the rectangles. As 
before by Whitehead's result, the Whitehead equivalence class of the resulting measured 
foliation is well-defined. 

If a measured foliation arises in this way from a measure on a train track, then we say 
that the train track carries the measured foliation. 

Let U(t) denote the cone of all measures on r, i.e., the subspace of K>q determined by 
the switch conditions. There is again the natural lR>o-action on U (r) — {0} by homothety, 
and the quotient V(r) = (U(r) — {0})/lR > o is the polygon of projective measures on r. 
Construction 14. II thus gives well-defined maps U(r) — > M.!Fq(F) and V{r) — > VJ-q{F). 

A recurrent train track is maximal if it is not a proper sub track of any recurrent train 
track. For general F = F£, complementary regions to a maximal train track are either 
trigons or once-punctured monogons, but in the special case of the once-punctured torus, 
a maximal train track has a single complementary once-punctured bigon. 

Theorem 4.2 jI5].|52*| For any maximal recurrent train track r in F , Construction ^ ■ 1\ 
determines continuous embeddings 

U(t) -> MFo(F) and V(t) -> VF (F) 

onto open sets. 

In fact, the geometric condition in the definition of train track precisely guarantees 
the injectivity in this theorem. In light of this result, one may regard a maximal train 
track in F as indexing a chart on the manifold VT§, and we next study the transition 
functions of this putative manifold structure. 

Torus Example For F = F\, two embeddings of train tracks as spine are illustrated in 
Figure 8, and in fact, every foliation is carried by one of these two train tracks. The cor- 
responding charts on the circle are also illustrated as well as the two points of intersection 
in the closures of these charts. 



Figure 8-charts for torus 

It is most convenient now to restrict to the "generic" case, where all switches of r 
are trivalent. For each switch of r, the decomposition of incident half-branches into 
incoming/outgoing thus consists of one singleton and one doubleton, and we say a branch 
of r is large if it is a singleton at both its endpoints (which are then necessarily distinct) 
as illustrated with the branch labeled e in the left-hand side of Figure 9. Likewise, if a 
branch is a doubleton at both its endpoints, then it is called small, while a branch which 
is neither small nor large is called half-large. 

Define the combinatorial splitting of a measured train track (r, /i) along a large branch 
e as illustrated in Figure 9, where we identify an edge with its /i transverse measure 
/i(e) for convenience. One imagines separating bands of horizontal leaves in the rectangle 
associated to e by excavating along the two "singular leaves" beginning at the endpoints 
of the large branch, i.e., beginning at the singular ties. If the measure fi is so that (either 
of) the singular leaves starting at an endpoint of e turn left or right, then the respective 
split is called a left (case 1) or right (case 3) split, while if the two singular leaves coincide 
for e, then the split is called a collision (case 2). 





b — a if a < b 



if a = b 



if a > b 



Figure 9-splitting 

Splitting and smooth isotopy of measured train tracks generates an equivalence relation 
on the set of all measured train tracks in F, and we shall let [r, /i] denote the equivalence 
class of the measured train track (r, //) and [r, p] denote the equivalence class of the 
projectively measured train track (r,p). 

Theorem 4.3 [32J //(r^/ij) are positively measured train tracks giving rise to corre- 
sponding measured foliations (J-i, fa) via Construction \4 ■ l\ fori = 1,2, then 

= [^2,^2] if and only if [n./ii] = [t 2 , fJ, 2 )- 

Thus, the space of all Whitehead equivalence classes of (projectivized) measured foliations 
is identified with the space of all splitting equivalence classes of (projectivized) measured 
train tracks up to isotopy. 

There is another aspect to the splitting equivalence relation on the set of all measured 
train tracks. In addition to splitting, one considers also shifting along a half-large branch 
by pushing two confluent branches of a train track past one another as illustrated in 
Figure 10. Shifting plays a role in the later discussion, and a basic result in train track 
theory [^2] is that if two train tracks are related by shifting, splitting, and smooth isotopy, 
then they are also related by splitting and smooth isotopy alone. 



Figure 10-shifting 

Let us finally give the idea of the proof that S is dense in MJ-q, as was mentioned 
before, by explaining density of S in each chart V(T) C "PJF , for some maximal train 
track r C F. We may approximate any fi £ V(r) C (M>o T ' ) — {0})/lR > o by a rational 
' ^ (Q>o — {0})/^>o (satisfying the switch conditions). Furthermore clearing 





measure ji 



denominators in //, there are N £ Z >0 and v £ (Z>[, r) - {0})/R >0 so that Nfi' 



v. 



We may construct an embedded family of curves in F from v by arranging u(b) > 
tie-transverse strands parallel to b in a tie neighborhood of r. By the switch conditions, 
there are at each vertex exactly as many incoming strands as outgoing, and there is a 
unique way to combine strands near vertices to produce a disjointly embedded family of 
curves. Let us give each component curve a weight 1/N and combine any parallel curves 
while adding their weights to produce our desired weighted family of disjointly embedded 
curves. Let [J 7 , /ix] denote the corresponding measured foliation (discussed before); tracing 
through the constructions, one finds that [r, //] = 1/N[t,u] gives rise to [J 7 , //i]. Letting 
\j! — > // and projectivizing, it follows easily that families of disjointly embedded curves 
are dense in V{r). With a little more work [32], one can approximate (in the topology of 
V(t)) such disjointly embedded families with a single curve, and this gives the asserted 
density of S itself. 

Thus, Thurston's boundary VJ r o{F) is a completion of the set S. In fact, one can 
approximate with a single non-separating curve (provided g 7^ 0); since any two such 
curves are equivalent under the action of MCg, it follows that the action of MC s g on 
VTq(F) has a dense orbit, and the maximal non-Hausdorffness of the quotient, which 
was mentioned before, is thereby established. 

4.3 Laminations 

Each basic formulation of the objects presented so far, namely, measured foliations and 
measured train tracks, requires passage to the quotient under an appropriate equivalence 
relation. Thurston has given a more ethereal, elemental, and elegant description of these 
objects as "measured geodesic laminations", where no passage to equivalence classes is 
necessary. Here we simply give the definition and a few basic properties referring the 
reader to [32] for instance for further details. A "lamination" £ in F is a foliation of 
a closed subset of F, and a "(transverse) measure" to C is defined much as before as a 
cr-additive measure on arcs transverse to L with the analogous condition of no-holonomy 
(where the homotopy is through arcs transverse to £ with endpoints disjoint from C). C 
is a "geodesic lamination" if its leaves are geodesic for some specified hyperbolic metric. 
The simplest case of a geodesic lamination is the geodesic representative of an element 
of S'(F). (In fact, for different choices of metric, the spaces of measured geodesic lam- 
inations are naturally identified via the circle at infinity in their universal covers, so we 
may speak of a geodesic lamination without the a priori specification of a metric.) A 
measured geodesic lamination has zero measure in F, and the intersection of a measured 
geodesic lamination £ with a transverse arc a in F is a Cantor set together with isolated 
points corresponding to intersections with simple geodesic curve components or arc com- 
ponents of £, if any. There is a natural topology on the set of all measured geodesic 
laminations in F, which is induced by the weak topology on the set of all 7i"i (F)-invariant 
measures supported on the Mobius band past infinity. M.C(F) (and M.Cq{F)) is the 
corresponding space of measured geodesic laminations (and with compact support) and 
corresponding projectivization VC(F) (and VCq(F)). A basic result in Thurston theory is 
MC(F) « MF(F), MC (F) « MF Q {F),V£{F) « VF(F), VC (F) « VF (F), where 
train tracks give suitable charts on any of these piecewise linear manifolds. Furthermore, 



the deformation theory due to Thurston, called "earthquaking" cf. [HHj, which we do not 
further discuss here, is most conveniently expressed in the context of laminations. 

4.4 Dynamics on train tracks 

In this section, we simply recall Thurston's classification of surface automorphisms as well 
as recall several basic facts about "pseudo-Anosov" mappings. In the process, we develop 
further basic techniques which will be required in quantization. Since Thurston's com- 
pactification produces a closed ball upon which the mapping class group acts continuously, 
one immediately is led to consider fixed points of this action. 

Theorem 4.4 Thurston's Classification ^3], [HE] Any orientation-preserving home- 
omorphism f : F — > F is homotopic to a diffeomorphism f : F — > F which satisfies one 
of the following conditions (and the only overlap is between 1. and 2.): 

1. /' fixes a unique point of T and is of finite order; 

2. /' is "reducible" in the sense that f fixes an element of S' ; 

3. /' is "pseudo-Anosov" in the sense that there is some A > 1 together with two 
measured foliations (J 7 ±,H±) , which share singular points and are otherwise trans- 
verse, so that f'(J 7 ±, jJ,±) = A ±x (J-'±,fJ>±)- The projective classes [F±,Ji±\ are the 
unique fixed points of f on T. The invariant A is called the dilatation of f or f. 

Notice the similarity with the trichotomy elliptic/parabolic/hyperbolic for fractional 
linear transformations. In the reducible case, one simplifies the dynamics by cutting F 
along a representative of the invariant element of S'. A pseudo-Anosov mapping is the 
analogue of an Anosov map of the torus in the current context of surfaces with negative 
Euler characteristic. 

In fact, train tracks provide a powerful tool for analyzing the dynamics of surface 
automorphisms owing to the fact that since a train track has a well-defined tangent line 
at each point, there is a coherent notion of a train track smoothly "carrying" a curve, 
another train track, or a lamination. 

Suppose that k is a smooth curve, a train track, or a measured geodesic lamination. 
We say that the train track r carries k and write k < r, if there is a C 1 map : F — > 
F homotopic to the identity, called the supporting map, so that </>(k) C t, where the 
restriction of the differential d<p p to the tangent line to k at p is non-zero for every p E k. 

We think of <fi as squashing together nearly parallel strands of k. For instance, any 
curve arising as before from an integral measure on r is carried by r, and more gener- 
ally, any curve, train track or lamination k which lies in a tie-neighborhood of r and is 
transverse to the ties satisfies k < r, where the supporting map collapses ties. 

For example, if a train track o arises from r by splitting and shifting (but no collaps- 
ing), then a < r, and we say that a arises from r by unzipping. 

Suppose that the train track o is transverse to a tie-neighborhood of r, say with 
supporting map : F — >• F . Let us enumerate the branches bj of r, for j = 1, . . . , n, and 
di of cr, for % = 1, ... ,171, and choose Xj G bj for each j. There is then an m x n- matrix 



A = (Aij) called the incidence matrix, where A^ is the cardinality of <f)~ 1 (xj) D cij. It is 
clear that the incidence matrix A : IR> — > IR™ describes the inclusion U{a) — > U(t) in 
the train track coordinates. 

We close this section with several basic results about pseudo-Anosov mappings. 

Theorem 4.5 A homeomorphism f : F —>■ F is a pseudo-Anosov map if 

and only if there is a train track t in F , with each component of F — r an at most 
once punctured polygon, so that r unzips to f(r) with no collisions. Furthermore, the 
incidence matrix A of the carrying f(r) < r is Perron-Frobenius, the eigenvector of A 
corresponding to the spectral radius A gives the projective measure fl + , and likewise the 
extreme eigenvector of the transpose of A gives 

Given a measured train track (r, //), consider the foliated neighborhood of r determined 
by via Construction 14.11 Choose some enumeration of the switches of r and serially 
follow the singular leaves from the switches until the first splitting (ignoring shifting), 
for the first switch, second switch, . . . , last switch, and then begin anew from the first 
switch. Suppose there are no collisions, and record the resulting sequence of right or left 
splits, so as to produce a semi-infinite word of rights and lefts. 

Theorem 4.6 [HH] The right-left sequence is eventually periodic if and only if the corre- 
sponding measured foliation is fixed by some pseudo-Anosov mapping. 

An explicit and simple construction of pseudo-Anosov maps is given by the following 
result. 

Theorem 4.7 [IT] Suppose thatC,D e S' admit representative arc families C,D inter- 
secting minimally which satisfy the condition that each component of F — U(C U D) is 
an at most once-punctured polygon. Take any composition w of Dehn twists to the right 
along elements of C and to the left along elements of T> so that the Dehn twist along each 
element of C orT> occurs at least once in w. Then w is pseudo-Anosov. 

Torus Example Consider a generic train track r in F^, so r has one large branch e and 
two small branches. The two small branches are canonically linearly ordered by first taking 
the branch a to the right and then the branch b to the left at either endpoint of e, and 
furthermore U(r) ~ M>q , i.e., the measures of the small branches a, b are unconstrained 

and uniquely determine the measure on e as well. Given a measure \i £ U(t) — {0}, start 
unzipping (r, fi) along either singular leaf, i.e., split along e, to produce another measured 
train track (ri,/xi); of course, t\ is combinatorially equivalent to r. For definiteness, 
suppose that B = fi(b) > fi(a) = A, so the split is a left split. The edge corresponding 
to b is the large edge of Ti, and the two small edges, in right/left order, have measures 
(A, B — A). Continue unzipping, i.e., next split (t\, [Ai) along its large edge to produce 
( T 2,l^2)- Again suppose that B > 2A for definiteness, so the second split is a left split as 
well, and the small branches of T2 have measures (A, B — 2 A). Continue unzipping (under 
the assumption that there are no collisions) until the first right split, say there are ai = 
[4] > 1 left splits before the first right split. Perform the right split along the large branch 



of ( 

T aii A*ai)j where the measures on the small branches are {A, B\) — (A,B — m\A\ to 
produce a train track whose small edges have measures (A — _B 1; Bi). Continue unzipping 
and suppose there are no collisions, i.e., suppose A and B are not rationally related, to 
produce a semi- infinite sequence of symbols L (for left splits) and R (for right splits). Let 
d\ denote the number of L's that begin this sequence, a 2 denote the length of the next 
consecutive sequence of i?'s, a 3 the length of the next consecutive sequence of L's, and so 
on. It follows from the discussion above that the continued fraction expansion of B/A is 
given by 



B/A = ai + 



a 2 + 



a 3 



Continued fractions occur in another related guise as well. The isomorphism MC\ ~ 

PSL-iffi) is induced by M 1— > jj and L h | ^ ^J, where M and L are the 

right Dehn twists on the meridian and longitude respectively. A mapping class in MC\ 
is pseudo-Anosov, reducible, periodic if and only if the corresponding fractional linear 
transformation is hyperbolic, parabolic, elliptic respectively. Every hyperbolic element of 
PSX 2 (Z) is conjugate to a product 

1 mA f 1 0\ (I m k \ ( 1 

1 y Vni 1 y/ "' \0 1 J \n k 1 

where rn^nj > are unique up to cyclic permutation. Furthermore, mi,ni, . . . ,mk,rik 
are the partial quotients of the periodic continued fraction expansion of the dilatation 
of the corresponding pseudo-Anosov map. It follows from this discussion that in MC\, 
all pseudo-Anosov mappings arise from the previous theorem. Indeed, the theorem gives 
the construction of two semi-groups corresponding to right/left or left/right twisting on 
meridian/longitude. For each semi-group, it is easy to construct a train track r in F± so 
that the matrix representation above precisely describes the action of the corresponding 
semi-group on the measures of the linearly ordered small branches of r; indeed, these two 
train tracks are illustrated in Figure 8. 



4.5 Decorated measured foliations and freeways 



In this section, we recall material from [33] which is required for quantization. If V C F 
is a cubic fat graph spine of F, then we may blow-up each vertex of T into a little trigon 
as illustrated in Figure 11. The resulting object r = t? has both a natural branched 
one-submanifold structure and a fattening, and furthermore, components of F — r are 
either little trigons or once-punctured nullgons. Thus, r is not a train track, but it is 
almost a train track, and is called the freeway associated to Y. Notice that each edge of Y 
gives rise to a corresponding large branch of r, and each vertex gives rise to three small 



branches. It is easy to see that every measured lamination of compact support in F is 
carried by the freeway r. The frontier of a once-punctured nullgon component of F — r 
is a puncture-parallel curve called a collar curve of F. A small branch is contained in 
exactly one collar curve, while a large branch may be contained in either one or two collar 
curves. 




Figure 11— freeway from fat graph 

A measure on a freeway r is a function /x 6 IR- 6 ^ satisfying the switch conditions, 
where we wish to emphasize that the measure is not necessarily nonnegative (as it is for 
train tracks). Let U (r) denote the vector space of all measures on r. Notice that /x G U(t) 
is uniquely determined by its values on the small branches alone, and the switch conditions 
are equivalent to the following "coupling equations" 

/i(ai) + /x(&i) = y.{e) = /u(a 2 ) + /x(£> 2 ), 

for any large branch e whose closure contains the switches V\ ^ t>2, where a^bi are the 
small branches incident on v,- t for % = 1,2. On the other hand, the values on the large 
branches alone also uniquely determine /x, and in fact, these values are unconstrained 
by the switch conditions. Indeed, letting cij denote the large branches incident on a little 
trigon with opposite small branches aj, for i = 1, 2, 3, we may uniquely solve for a measure 
/x on t, where 

K a i) = 7;{M a i) + M fl 2) + n(a 3 ) - %K a i)}i for * = 2 > 3 > 

and so we identify U(r) ~ M 1,5 ^), where LB{t) denotes the set of large branches of r. 

In particular, if \x is a nonnegative measure on r, then the analogue of Construction ^. ll 
in the current context produces a well defined equivalence class of measured foliations in F, 
where this measured foliation will typically contain a collection of puncture-parallel annuli 
foliated by curves homotopic to collar curves. Deleting these foliated annuli produces a 
well-defined (but possibly empty) class in MFq^F). Thus, a nonnegative measure on r 
canonically determines a point of AiF^F) together with a nonnegative "collar weight", 
i.e., the transverse measure of a transverse arc connecting the boundary components of 
the corresponding foliated annulus. 

In the general case that /x is not necessarily nonnegative, suppose that C is a collar 
curve of r. The switches of r decompose C into a collection of arcs, each of which inherits 
a corresponding real-valued weight from /x. Let {7i}i denote the collection of real numbers 
associated to the small branches of r that occur in C. Define the collar weight of C for /x 
to be fic — min{7j}™. Define a collar weight on F itself to be the assignment of such a 
weight to each puncture. 



We may modify the original measure \i 6 U(r) by defining //(&) = nip) — ji c if b is 
contained in the collar curve C for any small branch b of r. Thus, // is a nonnegative 
measure on /i with identically vanishing collar weights that determines a corresponding 
element of MFq(F). 

We are led to define the space M.Tq^F) = MFq^F) x IR s of decorated measured 
foliations and summarize the previous discussion: 

Theorem 4.8 [34J The space U{r) « M. LB ^ gives global coordinates on M.J-'q, and 
there is a canonical fiber bundle II : M.T§ — > M.T§, where the fiber over a point is the set 
IR S of all collar weights on F. 

Remark 4.1 The natural action of MC s g is by bundle isomorphisms of II. Furthermore, 
II admits a natural MC^-invariant section a : M.T§ — > M.T§ which is determined by the 

condition of identically vanishing collar weights. The restriction of a to AiFo C AiJ-'o 
gives a piecewise-linear embedding of the piecewise-linear manifold AiFo into the linear 
manifold (vector space) M.Tq « U{t) « IR 1 " 8 ^). 

4.6 Shear coordinates for measured foliations 

We now give an equivalent parametrization of measured foliations in terms of "Thurston's 
shear coordinates" that are close analogues of Thurston's shear coordinates Z a on Th(F). 
In fact, we have already encountered these quantities when describing the splitting pro- 
cedure train tracks (see Figure 9). There, excavating along two different singular leaves, 
we have obtained the "new" edge, which can turn either left or right (for splittings) or be 
absent (for collisions). 

We assign a corresponding signed quantity (positive for right, negative for left) as 
follows. Given a measure fi on the long branches of the freeway r associated to the 
fatgraph spine r C F, define the (Thurston's foliation-) shear coordinate of the edge 
indexed by a to to be 



in the notation of Figure 1 for nearby branches. ^From the very definition, ( a is inde- 
pendent of collar weights. Again, Thurston's foliation-shear coordinates are alternatively 
defined in terms of the signed transverse length of the arc between the singular leaves 
along Z a , in analogy to the geometric interpretation given before of the shear coordinates 
on T H . 

Note that the shear coordinates ( a are not independent. They are subject to the 
restrictions that 



for the sum over edges a & I surrounding any given boundary component, and we shall 
refer to these conditions as the face conditions for shear coordinates. Thus, the space of 
foliation-shear coordinates is of dimension LB(r) — n, where we let n denote the number 



C a = -^(A)-^B)+^C)-^D)) 




(4.1) 



of boundary components. One sees directly that for any assignment of shear coordinates, 
there is a well-defined point of M.T§ realizing them, thereby establishing a homeomor- 
phism between AiJ-'o and this sub-vector space M. LB ^~ n C IR LB ( T ) of shear coordinates 
on the long branches of r. 

To describe the action of the mapping class group on foliation-shear coordinates, we 
shall give the transformation under Whitehead moves, i.e., derive the analogue of formula 
()2.8j) for measured foliations, which is an elementary calculation using the formulas for 
splitting as follows. 

Lemma 4.9 Under the Whitehead move in Figure 2, the corresponding foliation- shear 
coordinates of the edges A, B, C , D, and Z situated as in Figure 2 are transformed 
according to formula \2.ty) 

M z : (Ca, Cb, Co, Cd, Cz) >-> 

(Ca + MCz), Cb ~ M-Cz), Cc + <Ph(Cz), Cd ~ 4>h{-Cz), -Cz) 



with 




Figure 12-foliation-shear coordinates 



Remark 4.2 Comparing expressions for the classical function <p{x) = log(l + e x ) and 
(|4.2jl . one finds that the latter is a projective limit of the former: 

<f) H {x) = lim \<p(\x) = lim ^ h (\x), (4.3) 

that is, all three transformations coincide asymptotically in the domain of large abso- 
lute values (or large eigenvalues for the corresponding operators) of Teichmuller space 
coordinates {Z a }. We shall actively use this property in Section 15. 3. 21 when proving the 
existence of the quantization of Thurston's boundary for the punctured torus. 



5 On quantizing Thurston theory 



5.1 Proper length of geodesies 

Definition 5.1 The proper length p. 1.(7) of a closed curve 7 in the classical or quantum 
case is constructed from the quantum ordered operator P 1 associated to a closed oriented 
edge-path with basepoint (to begin the linearly ordered word P 7 ) as 

p.l.( 7 ) = lim - tr log2T n (P 7 /2), (5.1) 

n^oo 77, 

where we take the principal branch of the logarithm and T n are Chebyshev's polynomials 
(cf. (|3.2()jl ). Since T n (cosh|) = coshy, it follows that p. 1.(7) agrees with half the 
hyperbolic length of 7 in the Poincare metric in the classical case. 

More explicitly in the operatorial case, we can determine p. 1.(7) explicitly in terms 
of the spectral expansion of the operator G 7 , which is known exactly. Namely, the basis 
of eigenf unctions of G 7 is "doubly reduced" in the sense that each eigenvalue (except 2, 
which is singular) with corresponding eigenfunction as, has the form e s ^ 2 + e _S//2 , where 
S ranges over the entire real axis, and as has the same eigenvalue as as. In fact, these 
functions coincide, so there is actually a representation on the positive real axis, which 
is nevertheless complete, and is singular at infinity and zero. We may define the proper 
length operator to the be one with the same eigenfunctions as for S positive (which 
constitute a basis in the function space) and with eigenvalues to be |5'/2|. This operator 
p. 1.(7) is then a well-defined operator on any compactum in function space. 

The proper length of a QMC or GMC C, again denoted p.l.(C), is the sum of the 
proper lengths of the constituent geodesic length operators (or the sum of half geodesic 
lengths calculated in the Poincare metric in the classical case) weighted by the number of 
appearances in the multiset. 

5.2 Approximating laminations and the main theorem 

Fix once and for all a spine r of F with corresponding freeway r. A measure /ionr gives 
rise to a (possibly empty) measured foliation in F together with a collar weight on the 
boundary components of F. Erasing collars yields an underlying measure \L\ > on r 
whose support is a sub-train track n C r, and the measured train track (n, fii) determines 
a (possibly empty) measured foliation. Via the canonical embedding of M.T§ into 
with vanishing collars, we may thus uniquely determine a point of AiJ-o by specifying 
foliation-shear coordinates on the long branches of r satisfying the face conditions 14.11 
i.e., AiJ-'o is naturally identified with a codimension n subspace of H LB ( T \ Passing to 
projective foliations, a point of VTq is given by the projectivization P( of a vector of 
foliation-shear coordinates £ = (Q), where £ G IR LjB(t ) — {0}, and i indexes the long 
branches of r, i.e., the edges of T. 

Definition 5.2 A sequence nP = (nf }, for (3 > 1, of integer-valued n.j, for i — 1, . . . , LB(r), 
on t is an approximating sequence for the projectivized measured foliation P( if the face 
conditions 14. II hold on nP and if lim^oo nf /n? = Q/Cj f° r a U h 3 with Q ^ 0. 



Constructed from n as an integral measure on r is a GMC C with integral collar 
weights. Just as with decorated measured foliations, components of C which are puncture- 
or boundary-parallel can be erased to produce a corresponding multicurve to be denoted 
Cfi- Ch is carried by a sub-train track of r, and it traverses the long branch of r indexed 
by i some number, say, rrii > of times, so rrii is the standard train track coordinate of 
integral transverse measure. In the usual notation as in Figure 12, one sees directly that 
nz = \{jtla — m B + m c - mo). We may also sometimes write for Ca 

Definition 5.3 A graph length function with respect to the spine Y is any linear function 

g\ s T {Cn) = g.\ S v {Crn) = W (5-2) 

i 

In particular, when all a« are unity, the graph length is just the combinatorial length of C^, 
i.e., the total number of edges of V traversed (with multiplicities) by all the component 
curves of C^. When the spine T and a are fixed or unimportant, then we shall write 
simply g.l.(Cy or g.l^C*^) for the graph length. 

Any graph length function is evidently additive over disjoint unions of multicurves, 
and more generally, is a linear function of m. 

We next describe the bordered or punctured torus case in detail. Each multicurve 
on the torus is uniquely determined by three nonnegative integers (mx,my,m z ) that 
satisfy one of the three triangle equalities rrii = rrij + m k , where {i,j,k} = {X,Y,Z}. 
Projectivization allows us to re-scale so that rrij and m k are relatively prime using this 
degree of freedom, so that the corresponding multicurve has just one component. 

As in Figure 8, the space VJ-q{FI) of projectivized measured foliations with compact 
support is a piecewise-linear circle for r + s = 1, and an alternative family of charts on this 
circle is given in Figure 13. The relation ~ in Figure 13 denotes the equivalence between 
different boundary cases between two different charts, and arrows represent one-simplices 
in VF{Fl). 
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Figure 13-the circle VC Q (Fl) 

In Figure 13, we use the previous notation Gx, etc. (see (|2.15jl ) but in a slightly 
different sense. Now, these quantities are (2 x 2)-matrices , not just geodesic functions, 
i.e., we do not evaluate traces in the corresponding formulas. There is thus an ambiguity 
in choosing the place in the graph where the matrix products begin. We indicate this 
place by drawing the reference cut (the dotted line). Changing the reference cut when 
moving along the circle corresponds to passing from one chart to another in the chart 
covering of the circle. Of course, choosing the reference cut does not affect the quantum 



trace operation. 

In order to have a good transition in the boundary cases, for instance, in the upper 
case in Figure 13 (and the other cases are similar and omitted), we must ensure that 
the corresponding functions for the quantities G^G^ and G™G% must coincide in the 
limit m — > oo with each other and with the corresponding quantity calculated merely 
for the "short" geodesic function Gx- To prove this, given two (2 x 2)-matrices Gx 
and Gz corresponding to geodesic curves, we can conjugate them by respective unitary 
transformations Ux and Uz to diagonal form with real eigenvalues e ±lx ^ 2 and e ±lz l 2 
since Gx,Gz are hyperbolic. We then have 



trG^=tr n V n V- 1 , V = U?U Z (5.3) 




and the proper length (j5.1|) is mix/ 2 + 0(1). 2 In order to have a well-defined projective 
limit, we shall "kill" the factor m in a consistent way, and this can be achieved by dividing 
the result by any graph length function of the curve since g\.(G^G z) = m ■ g.\.(Gx) + 
g-HGz). 

In the quantum case, however, the situation is much more involved. Indeed, let us 
consider an example of the product, which is of form U m V , as in (|5.3|h where U = e X//2 , 
V = e y / 2 , and [X,Y] = Anih. Thus, (U m )V = e ™x/2+Y/2+mirm/2 while ^([/™) = 
e mX/2+Y/2-mmh/2 ^ an( ^ ^ e corresponding logarithms do not coincide as m — > oo. More- 
over, even the Hermiticity condition does not often suffice to determine the proper length. 
For instance, given operators U = e aX and V — e tl39x , we may calculate that V{U m )V = 
e max +2i/3d x +maf3/3^ j e ^ proper length in this case is aX 2 + a/3/3 and depends on the 
parameter (3 (of course, this correction is purely quantum). This illustrates that proving 
continuity for the boundary transitions in the quantum case requires more subtle esti- 
mates, which we perform in the next section after deriving recurrence relations for the 
operators of quantum approximating multicurves. 

We may now formulate our main result on quantizing Thurston theory: 

Theorem 5.1 Fix a spine V of with corresponding freeway r. Fix any projectivized 
vector P( of foliation- shear coordinates on t and any graph length function g.l.. For any 
approximating sequence nP to P(, the limit 

lim P - L( ^ } (5.4) 
g.l.(C^) 

exists both in the classical case as a real number and in the quantum case as a weak 
operatorial limit. 

-i fa -r- 1 \ 

Unless the matrix U x GzUx has the form . For this matrix to determine a hyperbolic 

\r 4> ) 

element, the quantity <f> must be real greater than two. Multiplying this matrix by the diagonal matrix 

(0 — r" 1 e m ' x / 2 \ 

, so for sufficiently large m, the resulting product ceases 
re -rnl x /2 ( j )e -ml x /2 J 

to be hyperbolic, which is absurd. 



Because both the numerator and denominator in the limit are additive, this limit is 
projectively invariant and defines a continuous function (in the classical case) or a weakly 
continuous family of operators (in the quantum case) on the circle V£o(Fl r ), for r + s = 1. 

The proof of the previous theorem occupies the remainder of this section. The contin- 
ued fraction structure intrinsic to the torus case is used extensively, and various analogous 
operatorial recursions are derived and studied. There is a second essentially combinato- 
rial proof of this result, however, only in the classical case since we have no means to 
control the quantum ordering of the procedure. Nevertheless, the structures discovered 
are interesting, and we present this second proof in Appendix A, which depends upon 
the recursion (|5.8j) derived later in Lemma 15.21 Indeed, this basic recursion arises from 
"Rauzy-Veech-Zorich induction" [35] in the special case of the torus, which is derived 
from first principles in the next section. 

5.3 Elements of the proof 
5.3.1 Continued fraction expansion 

In each one-simplex in VCq(F), illustrated as arrows in Figure 13, the approximating 
multicurve is determined by two nonnegative integers, mi and m 2 , where we assume that 
rrii > ^2 with mi and m 2 relatively prime. It is convenient to represent the ratio m^jnii 
as a simple continued fraction: 

m 2 /mi = 

a i H 

a 2 + 

1 

+ r 

a n -i H 

a n 

We concentrate on the case nix = w*2> my — m ii m z = nix+niy, with the other cases 
following by symmetry, and describe the recurrence procedure for constructing the corre- 
sponding approximating multicurves (i.e., approximating geodesies, since any multicurve 
in the torus is just a multiple of a single geodesic). 

Referring to Figure 13, it is convenient to visualize this case by drawing a line 

© © 
• • • • o 

1 2 nil — 1 ^i ^i + 1 

which represents (the right side of) the cut over the edge Z (homotopic to the cut hz in 
Figure 12). The geodesic function is constructed as follows. Considering the multicurve 
in a neighborhood of the edge Z, we start from the lowest thread that goes from the 
edge Z to the edge X (labelled by the circled number one in Figure 12 and in (|5.6jl ) and 
moves to the right (this corresponds to the leftmost o in (|5.6jl ). We have the matrix Gy 
corresponding to passing consecutively through edges Z and X and come back to the line 



(5.5) 



o o 
nil + ^2 — 1 m i + ^2 



(5.6) 



(|5.6j) . which must be periodically continued, at the point that is situated to the right by 
m 2 points (in this case, just the leftmost • in (J5.6|) marked there by a circled number 
two). Each time passing the bullet sign, we must set the matrix Gx and passing the o 
sign we must set the matrix Gy. 

We thus move along the periodically continued line of circles and bullet signs in ()5.6|) . 
at each step jumping m 2 points to the right and setting the corresponding matrices Gx 
or Gy. We describe a recurrence procedure that produces from the continued fraction 
expansion ()5.5|) the correct sequence of operators, i.e., the correct ordered sequence of 
branches traversed by the corresponding curve. We shall split the construction procedure 
into stages. Let us consider the trajectory of the starting thread (the point m\ + 1 in 
(|5.6|) h Each stage terminates as soon as we come closer to the starting thread 1 than 
the distance to it at the beginning of the stage (approaching it from the opposite side 
compared to the beginning of the stage, see ()5.7|) below). We shall let Li denote the string 
of matrices for the first i stages, where by convention we arrange matrices from right to 
left. We terminate the first stage at the thread that is closer than m 2 to the left of the 
starting point, i.e., L 1 = (Gx) ai Gy. 

Define L Q = Gx and introduce Lj, for i > 1, which is the matrix Lj (composed 
from the elementary matrices Gx and Gy) in which the first two symbols of elementary 
matrices must be interchanged: if the first symbol of Li is Gy and the second is Gx, then 
the first symbol of Li is Gx and the second is Gy; all other elementary matrices retain 
their forms. We illustrate this procedure in Figure 14. 3 



H-l 



oooooooooooo 



end of ith stage 




end of (i — l)th stage 



end of (j + l)th stage 



(5.7) 



Figure 14— threads a curve 



3 Note that we terminate a stage whenever we come closer to the starting thread; this occurs at each 
stage on the opposite side from the previous stage, just as for continued fractions. The appearance of 
tilded quantities is explained as follow: each time during the recursion when we start a string "parallel" 
to some Lfe from the right to the starting thread (in the circled domain), we must follow the same string 
of branches because, by definition, there are no threads in the string except the very first thread that 
appears at a distance closer than \L^\ to the starting thread. On the other hand, if we start from the 
left of the starting thread, then we must interchange exactly the first two appearances of matrices in the 
resulting string. 



One thus sees directly that 



L 2 = (LO^LoLx 

l 3 = (1 2 y*- x l x l 2 

La — (L^) ai L 2 L 3 

l 5 = (uys- l LzU 



and this leads to the following recurrence relation. 

Lemma 5.2 Given the simple continued fraction expansion h5. <5|) of 1112/ mi, the sequence 
of matrices L n associated to the corresponding geodesic is given by the following recursion 

^ = (^)^" 1 ^-^ for i>l, L = G X , L\ — (Gx) ai G Y . (5.8) 
W = (L 2 i) 2i+1 L 2 i-iL 2 i 

Turning to the quantum case, we first show that the proper length operator (j5.1|) 
p.l.(*L 2i *)/g.l.(L 2 i) must agree with p.l.( * L 2i -x * ) / g.\.{L 2i -i) for a 2i -> 00, that is, the 
operators corresponding to continued fractions of form (|5.5|) with large coefficient a n 
must converge to the operator corresponding to the continued fraction terminated at 
the (n — l)st step. To this end, we must analyze the structure of matrix products and 
corresponding operators. Indeed, the operators Li and Lj enjoy elegant commutation 
relations as we shall next see. 

Notice that for every stage i we have a geodesic corresponding to the matrix Li (be- 
cause we can close the corresponding geodesic line without self-intersections), and we can 
therefore define the corresponding QMC 



d = * tr L 



't X 



The first observation pertains to Lj and Lf One of the corresponding curves can be 
obtained from the other by a parallel shift along the cut hz illustrated in Figure 12. 
Thus, the curves are disjoint and hence are homotopic on the torus, and so trLj = tr Li, 
i.e., 

C t = C u for i > 1. (5.9) 

Furthermore, the curves Li and can be perturbed to have exactly one intersection as 
one sees by considering how the corresponding geodesic curves pass through the cut h z .^ 
In the case where i = 2k, we obtain (in the notation of ()3.24|) ) 

[C 2k ,C 2k+1 ] q = HtiL 2k L 2k+l l, and [C 2k+1 , C 2k ] q = tr L 2k . 1 (L 2k ) a ^- 1 * , (5.10) 



4 One of these curves necessarily has odd subscript 2k + 1, and we can make a small shift of all threads 
of this curve to the right from the threads of the second curve L 2k (or L 2k +2 1 depending on the situation) . 
We find that there are no intersections of threads outside the region between terminating points of ^2fc+i 
and L 2k (or L 2k +2) 7 and in this domain, when closing the curves, there is produced exactly one intersection 
point. 



where the proper quantum ordering is assumed for the terms in the right-hand sides. 
Formulas (|5.1U|) are crucial when proving the continuity. Letting 



Im = x t r (-^2i-l) m 1 -^2i-2-^2i-l * , 

we find 

I m -iC2i-l = q 1/2 I m + q- 1/2 I m -2i (5.11) 

C 2 i-llrn-l = q- l/2 I m + q 1/2 Im-2i (5.12) 

where the basis of the recursion is given by Jo = e ly / 2 + e ~ ly l 2 and J_i = e' Z//2 + e ~ lz l 2 . 
^From (|5.11|) . we have the exact equalities 

I m = g- m /2j o(e m^/2 + g -mt x /2) _ g -m/2-l/2 g {m-X)l x /2 + e ~(m-l)i x /2^ (g ^ 

or, equivalently, 

J m = q m/2^ e me x /2 + e -m<x/2) /o _ g m/2+l/2 ( g ( m -l)^ /2 + g -( TO -l)^/ 2) (5U) 

We wish to present the expression (|5.13|) or f)5. 14j) in the form e mH ^+ H » ) where ifi, 
i^o are Hermitian operators independent of m, so the proper length is then just Hi (while 
H introduces quantum corrections that do not affect the proper length limit but are 
important for ensuring the proper commutation relations). 

Let us choose a compact domain T in the function space L 2 (IR) such that norms of 
all the operators in play are bounded for functions from this domain. Notice that neither 
term on the right-hand side in ()5.13|) is self-adjoint, but if one of these terms prevails 
(and it can be only the first term since we have an expression with coefficients that are 
all positive in the classical limit of the right-hand side, and hence we must have a positive 
left-hand side as well), then we can replace the total sum by this prevailing term and 
obtain an approximate equality 

g - m /2 /o(e ^x/2 + e - m W2) _ q m/2( e me x /2 + e ~ml x /2 )h 

in this limit. By considering the spectral expansion with respect to the eigenfunctions of 
the operator £x, we immediately conclude that I ~ e ~ 2mnd / d \ l x\ j n this limit, and then 

Hi = \£ x \/2, (5.15) 

where the modulus has to be understood in terms of the spectral expansion: having a 
QMC operator Cx which admits a spectral decomposition (see formulas ()5.24j) - (j5.26j) 
below) in functions \as), we define the operator \£x\ by its action on these functions: 

|**| |a s ) = |S| |a s ). 

A potentially problematic situation is when neither of the terms prevails and their 
difference remains finite as m — > oo. This would correspond, as in the classical case 
discussed before, to a situation where the corresponding element fails to be hyperbolic. 
As we next show, we must obtain "long" curves when m goes to infinity, so this is also 



impossible in the quantum case; that is, for Gx, Gy, Gz, and Gz from (|2.15|) . (|3.22|) and 
representing Gx = e £x ^ 2 + e~ ix ^ 2 , we must prove the operatorial inequality 

G Y e\ i *V 2 >q 1 ' 2 G z . (5.16) 

To prove this, first express e ix ^ 2 through Gx- Taking the positive branch of the 
square root, we find from the left-hand side of ()5.16|) the expression 



and since GyGx = q 1 ^ 2 Gz + q~ 1 ^ 2 Gz, we must compare two expressions Gy\J-f- — 1 

and (q 1 ^ 2 Gz — q~ l ^ 2 Gz) /2. Multiplying by the Hermitian conjugate on the right in both 
expressions, we eliminate the square root and arrive at Laurent polynomial expressions. 
After some simple algebra, we come to the inequality to be proved: 

-(qGzGz + q 1 GzGz) > G\. 

To see this, we expand the left-hand side 

{e -x-z + e x-z + e x+z + {q + q -i )e x ) 

+ q -±f^( e ~ X - Y - Z + e X+Y+Z +e z + 2e- z + (q + q' 1 )) 

+ q - t f^(e Y+Z + e ~ Y - Z + e ~ Y+Z + (q + q' 1 ) e - Y ) + ^±fl e ^ 

while the right-hand side is expressed as 

G 2 Y = e - x - z + e x ~ z + e x+z + 2 + (q + q~ l )(e x + e- z ). 

Subtracting this expression from the previous one, we obtain that this difference is 

q -^G 2 x + l±fL[e*+™ + e - X - Y - Z - q - rt (5.17) 

and both these terms are positive definite for \q\ = 1. 

The operatorial inequality f)5. 16j) has therefore been established. This proves that 
the limit (|5.4|) exists and is well defined at rational points of the continued fraction 
expansion, and we next turn to the case of infinite continued fraction expansions, i.e., 
infinite sequences of elementary operators. 



5.3.2 Mapping class group transformations and the unzipping procedure 

We consider now an infinite continued fraction expansion ai, a 2 , . . . , a n , a n+ i, . . . extending 
the notation of (|5.5jl . As we shall see, there is a corresponding sequence of unzippings 
of the freeway r associated to a spine r of F as in Section 4.5 as well as an associated 



sequence of mapping class group elements, expressed as Dehn twists, which reduce an 
approximating multicurve to one of two possible graph simple curves. 

Given the recursive representation ()5.8|) for the operator of a geodesic curve deter- 
mined by a continued fraction expansion ()5.5|) and applying two (unitary) operators Dx 
and Dy of the modular transformations of the form ()3.12|) that correspond to the respec- 
tive Dehn twists along the corresponding closed curves jx and jy (with the respective 
geodesic functions Gx and Gy), we shall construct the sequence of zipping or unzipping 
transformations. 

Definition 5.4 An approximating multicurve is determined by two nonnegative integers 
m\ and m 2 . As in Figure 13, enumerate such a pair as a triple (m 1 ,m 2 ,m 1 + m 2 ). If 
mi > m 2, the action of the Dehn twist Dy 1 along 7y (an unzipping transformation — the 
Dehn twist in the opposite direction) is 

Dy 1 : (mi, m 2 , mi + m 2 ) h-> (mi - m 2 , m 2 , mi) 

while if mi < m 2 , we apply the unzipping transformation along the curve 7x, which gives 

D x x : (nil, m 2i m i + m 2) l— * (^i; 1^2 — mi, m 2 ). 5 

Given a continued fraction expansion ()5.5|) . we construct the sequence of unzipping 
t r ansf or mat ions 

y (iorrnyori) ^y > ^o.io; 

which, when applied to the approximating multicurve (mi, m 2 , mi +m 2 ), reduces it either 
to (1, 0, 1) = 7y for n even or to (0, 1, 1) = 7x for n odd. 

Definition 5.5 Equivalently, we can consider the zipping procedure, that is, given a 
sequence of Dehn twists DpD^Dy 3 ■ ■ ■ D^' 1 D a x applied to the curve j Y , we obtain the 
curve (mi, m 2 , m x + m 2 ). 

Considering the sequence (j5.18J) of quantum Dehn twist operators ()3.12|) and exploiting 
the quantum invariance from Lemma f3. 91 we come to the main observation that having 
an involved expression for the proper limit (|5.1|) of a QMC operator constructed by the 
rules described in Lemmas 13.91 and 15.21 in terms of the elementary operators X, Y, Z, we 
may perform the sequence (j5.18|) of unzipping quantum modular transformations, which 
reduces this operator to a standard form of the quantum operator Gy or Gx expressed 
through the new operators X^ n \ Y^ n \ related to the initial operators by this sequence 
of quantum modular transformations. This is the operatorial statement of naturality of 
lengths under the mapping class group action. 

It is intuitively natural to imagine that as the geodesic lengths must diverge as mi, m 2 
tend to infinity, we must eventually come to an asymptotic regime where all quantities 
are large in the literal or operatorial sense for all sufficiently large n. The 
inexorability of the approach to this asymptotic regime is not obvious and is described in 
the next section. 

5 In terms of the symbolic dynamics of elementary operators Gx and Gy , we can present this action 
as follows: D^(G X ) = G x , D x x (G x Gy) = Gy and D Y 1 {G X G Y ) = G x , D Y 1 {G Y ) = Gy. 



5.3.3 Asymptotic regime 



Let us recall the modular transformations for X, Y, and Z variables: 

D~ l (X, Y, Z)^(X + 2<j) n {Z), -Z, Y - 2<j) h {-Z)) (5.19) 

and 

Dy 1 (X, Y, Z) ^ {-Z, Y - 2<j) h (-Z),X + 2<\>\Z)). (5.20) 

In terms of the quantities U = e X//2 and V = e ~ Y ^ 2 in the case where X + Y + Z = 0, 
we have 

e x/2 + e -x/2-y \ fu + VU- 1 v\ , , 

°i I = I e -v/*^ H ^tr'v ' (5 ' 21) 
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Df UJ y = U-^+^J s U+^J- (5 - 22) 

Worth mentioning is that since each operator Z a is Hermitian, each exponential is positive 
definite, so we can always write, for instance, that U +VU~ 1 V > U in the sense of spectral 
expansion: (f\U + VU^Vlf) > (f\U\f) for any function / G L 2 (R). 

Using now an alternating sequence of transformations (j5.21j) . (|5.22j) . we shall subse- 
quently show that we attain the asymptotic regime of large positive X (large U) and large 
in absolute value negative Y (large V) starting from every pair of X and Y lying in a 
compact domain of the (X, F)-plane in the classical case or acting within a compactum of 
test functions with bounded derivatives in the function space in the quantum operatorial 
case. 

In this section, we verify that the asymptotic regime is attained for the distinguished 
sequence of modular transformations corresponding to the Fibonacci number sequence 
(golden mean), namely, for alternating D x l and Dy 1 . The proof in the general case is 
analogous although the structure is more involved, as described in the next section. 

Given the sequence of transformations D^Dy 1 • ■ ■ D^Dy 1 =(D] c 1 Dy 1 ) n , we obtain 

Dx'Dy 1 ( U y J = (5.23) 

UV^U+V 

V^tV^Uq- 1 ' 2 + U^Vq^UiV^Uq- 1 / 2 + U^Vq 1 / 2 ^ 1 / 2 + V^U^V 1 ' 2 

In the classical case, the asymptotics is already clear from this formula; for the first 
entry in (|5"23Jl . we have UV^U + V = UiV^U + U^V) > 2U as the expression in the 
parentheses has the form e 5 + e~ s > 2 for any real S. The same logic applies to the 
second entry in ([5.23)1 . and we deduce that the classical part V(U + U' 1 ) has the same 
property. 

The proof given below in the quantum case is more subtle as it needs a thorough 
operatorial analysis. Nevertheless, the estimates turn out to be close to those in the 
classical case, which we briefly discuss here: we must prove that a lower bound on the 



operatorial spectrum on a compactum in the function space diverges with n. This is a 
routine procedure, which uses that the action of operators X and Y in the basis of, say, 
normalized Hermitian functions h n has the form 

X\h n ) = \f<bxh(\/n + l\h n+l ) + y/n\h n -i)) 
Y\h n ) = V47rh^(Vn+ l\h n+1 ) - y/n\h n ^i)). 

These operators "almost" commute in the domain of large n, which allows the combina- 
torics to be analyzed semiclassically. 

In the quantum case, we recall the construction of quantum Dehn twists and their 
eigenfunctions from [2E]- The generator of the Dehn twist Dx has the form 

D x = e^ 2mh F h (q 1+Pl ), Ql = X/2, Pl = 2mhd x (5.24) 

and because it commutes with the geodesic length operator Gx, they share the common 
set of eigenfunctions 

\a s ) = e ~^ 2 n^F\S + X)F h {-S + X) (5.25) 

with the eigenvalues 

G x \a s ) = 2cosh(5/2)|a 5 ); D x \a s ) = e s2 ^ ih \a s ). (5.26) 
The functions \as) constitute a complete set of functions in the sense that 
(a T \a s ) = 5{S - T)v-\S), u(S) = 4 sinh(TrS) smh{nhS) 

and 

poo 

/ v(S)dS \a>s) (as\ = Id 
Jo 

We now split the plane of the variables (X,Y) = (X^,Y^°>) into four sub-domains 
and consider the action of the Dehn twists D x l and D Y l in each sub-domain. 

Domain I. {X<°) > 0, Y^ > 0} U {X<® > 0, < and |X<°)| > |F (0) |}. 
D x l (X(°),F(°)) = (X«y«) = pf(°),Jf(°> +y<°>) e Domain I; 
D Y l (X(°),y(°)) = (X«,y«) = (X^ +yW,-2X(°) -yW) e Domain II; 

Domain II. {X<® > 0, Y<® < and \X<®\ < |y (0) |} U {X<® < 0, Y<® < 0}. 
D x x (X(°),y(°)) = (X« y«) = - 2y(°),X(°) + Y®) e Domain I; 

/V (X(°),y(°)) = (XW,y«) = pr<°> +y< >,y(°>) e Domain II; 

Domain Ilia. {X^ < 0, y<°) > and |XW| < \Y (0) \}. 
D x l (Jf(°), yW) = (X«, y«) = (XW, X^ + yW) e Domain Ilia or Illb; 
D Y l (X(°), yW) = (X«, y«) = (XW + yW, -2X(°) - yW) e Domain I or II; 

Domain Illb. {I<°) < 0, Y^ > and \X<® \ > \Y^\). 



D x l (X< >,Y<°>) = (X<V,Y<U) = -2Y®,X® +y(o)) G Domain i or II; 

Dy 1 (XW,^^) = (lW,yW) = (XW +y(°),F( )) G Domain Ilia or Illb. 



We see that only Domain III is potentially problematic. This regime is however un- 
stable: absolute values of X and Y variables decrease in this regime and they eventually 
leave the asymptotic regime as soon as we remain in Domain III; immediately upon leav- 
ing this domain, we come to domains I and II and will never leave this three quarters 
of the (X, K)-plane. The above considerations of U, V just demonstrate that even if we 
were initially in Domain III, we come to the nonasymptotic domain of bounded X and Y 
and then will leave this compactum moving toward asymptotic expansions in domains I 
and II. 

The asymptotic dynamics always takes place in the first three quarters of the (X, Y)- 
plane. Nevertheless, even this dynamics is rather involved. The stable regime corresponds 
to the case where we are in Domain II before applying one or several operators D x l . The 
application of the first of these operators brings us to Domain I, and upon subsequent 
applications of the operators D x x we remain in Domain I. Next, if we were in Domain I, 
then the very first application of the operator Dy 1 brings us to Domain II, and we then 
remain in Domain II upon subsequent applications of Dy 1 . 

We turn now to actual geodesic lengths of curves or proper lengths of operators. If 
a sequence of unzipping transformations terminates, this means that we have a graph 
simple geodesic, which is either Gx if the last transformation was Dy 1 or Gy if the last 
transformation was D^ 1 . Considering the corresponding geodesic or proper lengths, we 
find that up to exponentially small corrections, the leading contributions in the above 
domains are 

p.l.( 7y ) = X + Y/2 in domain I, (5.27) 
p.l.( 7x ) = -Y-X/2 in domain II (5.28) 

(see expressions (12.15(0 . 

Thus, although the transformation laws for the variables X, Y themselves do not 
possess the property of linearity with respect to the parameters <Zj, a,j, when applying 
sequences of transformations (D~ x 1 )1 = D x ai or (Dy 1 )^ = D y a \ the proper lengths do 
possess this property! Namely, starting with variables (X^°\ Y(°>) lying in the correspond- 
ing domains and applying the sequences of transformations D x ai or D y a \ we obtain for 
the resulting proper lengths the following expressions: 

y(o) / ls v(o)\ 
p.l.( 7yw ) = _ + a t ( -F (0) — J for X(°), G domain II (5.29) 

v(o) / , y(o) \ 
p.l.( 7x0) ) = — + aj ( X (0) + — J for X(°), G domain I (5.30) 

Let us now explore the asymptotic formulas ([5 .27)1 . ()5.28|) and (j5.29j) . ()5.30|) first in 
the classical case to close this section, relegating the discussion of the quantum case to 
the next section. 



Assume that we start from the variables (X^°\Y^) in domain I and have the cor- 
responding initial length p. 1. (7y(o)) from (|5.27jl . Applying the transformation D y a \ we 
obtain new variables (XU),Y®) and the new proper length P-L(7xW)) ()5.30j) having form 
f!5.28|) in these new variables, which must now lie in domain II. Note that explicitly 

{X®, Y®) = (X® + F (0) + (aj - 1)(-2X<°> - F (0) ), - F (0) ). 

We then apply the transformation D x a% to obtain variables (X^' 1 ', Y^'^), and the proper 
length p.l.(7y(j,i)) is expressed as in (J5.29|) . where the term multiplied by is none other 
than p.l.(7 X (j)) and the term — is exactly p.l. (7y(o)). We thus find in the asymptotic 
regime that the corresponding lengths are related by exactly the same recurrence relation 
as for a graph length (the latter of which follows immediately from (|5.8jl ): 

p.l.(7yCi,0) = Oip.l.(7 x U)) +p.l.(7y(o)), (5.31) 

g.l.(7 Y y,i)) = aig.\.{j x u)) +g.l.(7y(o)). (5.32) 

It is then easy to conclude that the ratio of these two quantities has a definite limit 
as i — > 00 for any sequence of numbers aj. It is a standard estimate: given two numerical 
sequences (J5.31|) and ()5.32j) and denoting the relative error of their ratio cLS 1. e., at 
the ith step, the ratio is S(l + £j), where S is constant, for a^+i > 1 at the {i + l)th 
step, we obtain < £j/(aj+i — 1/2), or if we have two coefficients a i+ i = a« = 1, then 
Ei + \ < £j_i/1.5. In general, for e, small enough, we always have e.- l+ i < Si. This shows 
that the relative error decreases exponentially with the index i. 



5.3.4 Quantum continued fraction expansion 

Let us turn again to the sequence (|5.18|) of unzipping transformations. In order to obtain 
operatorial expressions, we consider the unitary operators D x , Dy and explicitly indicate 
the variables in which these operators are expressed, i.e., we write D xU) =D X (X^,Y^) 
for the Dehn twist along 7x at the (j + l)th step. Thus, 

fyW y{n)\ _ r)-a n n-«»-i ... n-«2 n- fl ify(0) y (Oh r^i n«2 . . . n"n-i rta n 

(5.33) 

In order to represent such long strings of operators in terms of the original operators 
(X^°\Y^), we invert the dependence, i.e., we remember that, for instance, 

7")— 0,2 7-)— ai 7-)— C12 r-jai 

xf 1 ) ~ x(°) YW' 

etc., which gives 

frW — n~ ai n~ a2 ■ ■ . n-"-™- 1 n- Q n (x^ v ( -°hn an n an - 1 . . . n a2 n ai ia\ 

\^ > 1 1 — ^yw x(°) y(o) x(°)'V ' / x(°) y(°) x(°) y(°) V u - J ^7 

We shall compute with bases of functions that are convenient in the asymptotic regime. 

Let 

!/„,,) = e^" s ) 2 / 2 fx,sER. (5.35) 
These functions constitute a basis at each /x: 

(U,t\U,s) = S(s-t), / ds\f lh .){f„\ = — Id. (5.36) 



For two arbitrary real numbers w and 7, we have 

e iwx2/2 \U,s) = e^ 2 / 2 f^,^)- (0.37) 



, ildi/2 



V 1 + 7/^ 



!/„,.) = _ / 1 V (5.38 



and 



(/u>|/u,t> = —J ! {s-t) + —s5{s-t) ) (5.39) 



z/i z /i 

(./;,, ./;,/) = — (5.40) 



Z Z/i 

We now define the dimensionless variable x and set 



-1 d 



X\f^) = VAtcHx ■ 1/^); F|/ M)S ) = V4^-^ • |/ M , S ). (5.41) 

The explicit formulas for the operators D x and D Y acting on \f^ s ) in the asymptotic 
regime are 

D^f^s) = e^-V* 2 ' 2 e~ id * e*^/*.). (5 . 42 ) 

£>y%, s > = e^- 1 ^ 2 e"» 2 e ia2 / 2 |/ M , s ). (5.43) 

In order to establish the required recurrence relation, we must compare matrix elements of 
the three consecutive length operators in the corresponding operatorial decompositions: 

4? = (u, s \x+^\f, >t ), 

4? = (U,sD Y a >\-Y-j \Dpf,, t ), 
4f = (U, s D Y a >D x a >\X + ^\D^f,, t ). 
Now, using formulas (|5.37|) - (|5.43j) . it is straightforward to show that 

A<* t > = AV>+a i A<£ , (5.44) 

for all s, t, i.e., we again attain the recurrence relation ()5.3H) in the asymptotic regime, but 
now for the matrix elements of the operators of the quantum proper lengths. Estimates 
show that the corrections due to both the (operatorial) deviations from the asymptotic 
regime and the error parameters £j (as for (j5.31|) . (j5.32|0 decrease exponentially with 
the index z, so the limit ()5.4|) exists in a weak operatorial sense. We conclude that 
ratios ()5.4|) define a weakly continuous family of operators parameterized by projective 
transverse measures on the freeway associated to a spine of the once-punctured torus. 
This completes the proof of Theorem 15.11 



6 Conclusion 



We hope to have added to the mathematical foundation and general understanding of 
the quantization of Teichmuller space and its geometric underpinnings in the first several 
sections of this paper. We also hope that the survey given here of train tracks and their 
extensions might be useful. 

The quantization of Thurston's boundary in general seems to be a substantial project, 
which we have only just begun here with the quantization of continued fractions. First 
of all, one would like a better understanding of the operators we have constructed, for 
instance, an intrinsic characterization or an explicit calculational framework for them. At 
the same time, our current constructions depend upon a choice of spine, and there would 
seem to be a more invariant version of the theory, where the choice of spine is dictated by 
the combinatorics of the cell decomposition of Teichmuller space; the calculations in this 
paper apply to each such spine (since there is a combinatorially unique cubic one) for the 
once-punctured torus. 

Second of all, the quantization of Thurston's boundary for higher-genus or multiply- 
punctured surfaces may be approachable using the improved quantum ordering. Namely, 
in any fixed spine of the surface, there is a fixed finite family of "edge-simple" closed edge- 
paths which by definition never twice traverse the same oriented edge. It is elementary 
to see that any closed edge-path on T may be written non-uniquely as a concatenation 
of edge-simple paths, where the particular concatenation depends upon a starting point. 
(Edge-simple paths were studied as "canonical curves" on train tracks in |8.; they contain 
the extreme points of the polyhedron of projective measures on the track.) It follows 
that an arbitrary leaf of a measured foliation carried by a freeway can be written as a 
concatenation of paths from this finite collection of edge-simple paths. The corresponding 
quantum operatorial statement results from the improved quantum ordering described 
here. Thus, whereas the quantization of the once-punctured torus devolved, in effect, to 
an analysis of two-letter words, the quantization of Thurston's boundary sphere in general 
may depend upon an analysis of words comprised of letters which are edge-simple paths. 

One appealing long-term goal would be to discover the Thurston classification already 
on the operatorial level, for instance, with the dilatation in the pseudo-Anosov case ex- 
plicitly computable from the MCG operator or from the invariant projective foliation 
operator. 

Another intriguing aspect involves generalizations of graph length functions insofar as 
the proof of Theorem 15.11 holds taking as graph length any continuous positive definite 
function which is homogeneous of degree one. A natural choice of such a function is 
induced by the geodesic length of the corresponding geodesic curve taken for a fixed 
basepoint in Teichmuller space on its fixed spine, for instance, vanishing shear coordinates 
on the usual spine in the once-punctured torus. What sort of regularity (e.g., piecewise 
smoothness) is achieved in the operators corresponding to points of Thurston's boundary 
under such "gauge fixing"? 

Also worth mentioning are very recent advances in the description of quantum sl(n, K) 
connections [40 J, where one finds an improved quantum ordering in a more complicated 
higher- dimensional setting. 



Appendix A Combinatorial proof of Theorem 15.11 

In this appendix, we give a complementary, combinatorial proof of the classical Theo- 
rem |0] using the recurrence relation (|5.8j) . At the present state of understanding, the 
proof applies only to the classical case as we cannot control the quantum ordering. 

Let us recall the structure of the matrix product (|2.5jl . It is a sequence of matrices 
Lz, Rz with different Z. It can be always segregated into clusters of matrices 

Lg = L Zl+s L Zt+s _ 1 ■■•L Zi 

and 

Rz = Rz,.J-z,., , ■ • • L Zy 

The periodic extension of expression ()5.8|) is always an alternating sequence of matrices 
Lf? and Rg\ 

Pzx,...z n = ■■■ {Lz a Rz s ^) ' ■ ' ( L z R z.) 

First note that it is impossible to have arbitrarily long sequences of only left or right 
matrices for a given graph: the maximum length is restricted to be less of equal the 
maximum graph length of geodesies around holes. This means that the length of a single 
cluster for a given graph is always bounded once the topology is fixed. 

One can directly calculate the product (Lg Rg ) for Lg = L Zl ■ ■ ■ L Zm and Rg = 

Rz m+ i ' * ' Rz m+k - 

{LgRg 2 ) = As+s+ + + SisJ + S 1 S 2 ) 

+D(s+si +stS 2 ) + PS lS +, (A.l) 

where Sj, sj, for j = 1,2, are the following coefficient functions: 

s f = e ±Er =1 V2 ; a ± = e ±sr= + m fe +1 v2 

m m+k 

S 1 = e + ^ Zr/2 ~^ Zr/2 , S 2 = e-^r=^ z ^ 2+ ^ q k z r /2^ 

q=2 q=m+2 

and A,B,D,P are the special (2 x 2)-matrices ("letters"): 

B = 





P= . „ ■ (A.2) 



These letters possess interesting multiplication properties which are summarized in the 
next lemma, whose proof is a routine calculation. 



Lemma A.l The alphabet lemma. The multiplication table of letters (jA.2|) reads: 





A 


B 


D 


P 


A 


A 





D 





B 


P 


B 


B 


P 


D 


A 


D 


D 


A 


P 


P 





B 






(A.3) 



so the trace of any product of these matrices is either unity or zero. In the product of t 
matrices of form (1A. ljl . the only monomials that survive are 



[{A + D) ia D]B j e(B + P) ■ ■ ■ [(A + D) lp D}B^(B + P) and {A + D)\ E 1 



(A.4) 



The main point is that almost all cancellations of letters in long words are due to the 
local multiplication rules (|A.3J) . This means that, having a long sequence of letters, say, 
Li + n from Lemma E21 we can split it into pieces depending on sequences of letters Li, 
L/_i, Li, and L/_i, where the index / is also assumed to be big enough. That is, let 
Li+n = LiLi-iLiLi . . . Li comprise pn entries Li and Li and qN entries L/-i and Lj_i. 
We have then the following estimate: 6 



log tr Li +N - p N log tr L/ - q N log tr Lj_i | < C (q N + p N ) , 



(A.5) 



where the constant C depends only on the Teichmiiller space coordinates Z a and on the 
genus and the number of holes of the Riemann surface, and we have also used (J5.9|) . The 
ratio of the coefficients is given by the continued fraction 



Qn/p 



1 



N 



+ 



1 



(A.6) 



+ 



1 



o>i+n 



and also has a definite limit as — ► oo. Now the estimate follows: up to exponential 
corrections, p.l.(L) coincides with the log trL, so for any e > 0, let us choose the index / 
such that ep.l.(L/_i)/2 > C and ep.l.(L/)/2 > C. Thus, 



p.l.(L)(W; 



I+Nj 



:i + 0(e/2)) 



pjyp.l. (L)(Li) + q N ^.\.{L){Li^) 
p N g.\.(Li) + givg.l.(L/_i) 



(A.7) 



and because the ratio qN/PN has a definite limit as A^ — > oo, there exists N such that 
the relative error of this ratio times the sum of ratios of proper and graph lengths of L/ 
and Li_i will not exceed e/2. Thus, the collective relative error for such fixed / and for 
all A^ > Ao is less than e, proving the theorem. 

6 This estimate also follows from the properties of long geodesic lines in hyperbolic geometry: for two 
lines of large lengths Li and L 2 intersecting at angle a, the length L 3 of the third side of the resulting 
triangle is L\ + L 2 + log((l — cosa)/2) + 0(1/L). This also shows that our estimate is very rough. 



Appendix B Degeneracy of the Poisson structure 



We shall explicitly calculate the degeneracy of the Poisson brackets ()2.9|) for a special 
graph and choose the graph whose "building blocks" are depicted in Figure 15. Namely, 
we have a line tree subgraph comprising edges Xj with attached subgraphs as in Figure 15a 
and 15b. Attaching a subgraph of type a corresponds to adding a handle (increasing g by 
unity) while a subgraph of type b corresponds to adding a hole (increasing s by unity). 
We shall assume that 2g + 2s > 5 to avoid the once-punctured torus, which is already 
handled separately in Section 3. 



A Ei Nj 




a b 
Figure 15-building blocks 

For the variables A, B, C, D, E in Figure 15a, we have the Poisson bracket (sub)matrix 





A 


A 


Q 


A 


Ei 


A 





1 


-1 








Bt 


-1 





1 


1 


-1 


Ci 


1 


-1 





1 


-1 


A 





-1 


-1 





2 


Ei 





1 


1 


-2 






where the entries are the Poisson brackets between the corresponding variables. Adding 
the last row to the next-to-the-last row as well as adding the last column to the next-to- 
the-last column, then adding the third row to the second row as well as the third column 
to the second column, we obtain the matrix 



( 





-1 





o\ 








1 





-2 


1 


-1 








-1 














2 


v° 


2 


1 


-2 





which obviously has rank four and can be further reduced (without adding the first column 
or row to any other) to the form 



/ 

















\ 










+1 















-1 



























+2 




V 











-2 





/ 



Thus, erasing all columns and rows corresponding to the variables P>i, Ci, Di, and Ei 
leaves invariant the rank of the Poisson bracket matrix. 

Adjoining the subgraph in Figure 15b creates exactly one degeneracy as the variable 
Nj Poisson commutes with everything (as it must be when adding a hole). 

It remains only to calculate the rank of the matrix corresponding to a tree graph with 
edges Xi and A, remaining after erasing all B-, C-, D-, E-, and iV-variable rows and 
columns. The corresponding Poisson bracket matrix has dimension 2g + 2s — 5 > and 
the simple block-diagonal form 






1 -1 









-1 
1 


1 
-1 




1 -1 








-1 
1 


1 
-1 




1 -1 








-1 
1 


1 
-1 















Adding each even-index row to its predecessor as well as adding each even-index column 
to its predecessor, this reduces to the matrix whose only nonzero elements are +1 on the 
main super-diagonal and -1 on the main sub-diagonal. Since this matrix has full rank, 
the discussion is complete. 
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